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Abstract 

The universal Kummer threefold is a 9-dimensional variety that represents the total 
space of the 6-dimensional family of Kummer threefolds in P 7 . We compute defining 
polynomials for three versions of this family, over the Satake hypersurface, over the 
Gopel variety, and over the reflection representation of type E7. We develop classical 
themes such as theta functions and Coble's quartic hypersurface using current tools 
from combinatorics, geometry, and commutative algebra. Symbolic and numerical 
computations for genus 3 moduli spaces appear alongside toric and tropical methods. 
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1 Introduction 

Kummer varieties are quotients of abelian varieties by their involution x 1— > —x. Each g- 
dimensional Kummer variety has a natural embedding into P 29 " 1 by second order theta 
functions. The moduli space is recorded in a second copy of P 29 " 1 by the corresponding 
theta constants. By taking the closure, this construction defines the universal Kummer 
variety K g . This is an irreducible projective variety of dimension ( 9 ^ 1 ) + g, defined over the 
integers, and naturally embedded in P 29-1 x P 29_1 . The image of the projection of JC g onto 
the P 29 " 1 of theta constants is the Satake compactification of the ( 9 1 2 1 ) -dimensional moduli 
space A g (2,4). The various Kummer varieties of dimension g appear as fibers of this map. 

Our object of interest is the universal Kummer ideal X g . This is the bihomogeneous prime 
ideal in a polynomial ring in 2 9+1 unknowns that defines JC g as a subvariety of P 29 " 1 x P 29 " 1 . 
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What motivated this project was our desire to understand the ideal X 3 of the universal 
Kummer threefold in P 7 x P 7 . Grushevsky and Salvati Manni write in [28, §6] that the ideal 
X 3 "is known quite explicitly", and we were wondering how to communicate its generators 
to Macaulay 2 [27]. Our current state of knowledge about this question is presented in 
Section 8 of this paper. The example of X2 is worked out in Example 1.1 below. 

In our attempts to understand X 3 , we also studied a variant of /C 3 , over a base Q that 
is a moduli space for plane quartics. We discovered a beautiful mathematical story that 
connects classical algebraic geometry topics, such as Coble's quartic hypersurface, Gopel 
functions, and the type E 7 reflection arrangement, with more modern topics, such as toric 
geometry, tropical methods, and numerical computation. Our study of this variation was 
originally inspired by Vinberg's theory of representations [52], but can be defined without 
any reference to this theory. However, as we shall see in Section 8, this point of view gives the 
construction of a large number of equations that we cannot see how to construct otherwise. 

Moduli of plane quartics have been studied extensively, notably in [11] and [18, §IX]. A 
key feature of Kummer varieties of the Jacobians of plane quartics is this: in their embedding 
111 P 7 , there is a unique quartic hypersurface, called the Coble quartic, whose singular locus is 
the Kummer variety. We are interested in the moduli space of Coble quartics. This moduli 
space can be parametrized by the Gopel functions mentioned in [18, §IX]. These functions 
embed it as a subvariety Q of P 134 , so we call it the Gopel variety. In fact, it sits in a 
linear P 14 C P 134 , and it can be alternatively parametrized by a Macdonald representation 
of the Weyl group of type E7. Hence the rich combinatorics of reflection arrangements is 
embodied in Q. Colombo, van Geemen and Looijenga [12] studied the Gopel variety from 
this perspective, and realized it as a moduli space of marked del Pezzo surfaces of degree 2. 

It is natural to ask for the prime ideal of Q embedded in P 14 . We solve this problem and 
calculate its graded Betti numbers. A delightful picture emerges in the embedding in P 134 . 
Together with trinomial linear equations that define P 14 in P 134 , all equations are written as 
binomials. This shows that there is a larger toric variety T sitting in P 134 that contains Q 
as a linear section. In fact, there are 35 cubics and 35 quartics that cut out Q. The cubics 
are genus 3 analogues of the toric Segre cubic relations of Howard et al. in [32, (1.2)]. 

We obtain a range of results on T and the corresponding 35-dimensional polytope A. It 
has 135 vertices and 63 distinguished facets, and their incidence relations admit a compact 
description in terms of the finite symplectic space (F2) 6 . Furthermore, it exhibits symmetry 
under the Weyl group of type E 7 . It should be very fruitful to further study the combinatorics 
of the polytope A, such as its face lattice and Ehrhart polynomial. Even more important, 
this also paves the way for studying the tropicalization of Q, which we hope will provide a 
useful model for genus 3 curves over fields with a non-trivial, non-archimedean valuation. 

All relevant basics regarding theta functions and abelian varieties will be reviewed in the 
next two sections, along with pointers to the literature and to software. As a warm-up, we 
first present the solution to our motivating problem for the much easier case of genus g = 2. 

Example 1.1 (The Universal Kummer surface). The five- dimensional variety K.2 is a hy- 
persurface of degree (12,4) in P 3 x P 3 . Its principal prime ideal X 2 in the polynomial ring 
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Q[uoO) w oi> u w, u u] x oo, x oi, x io, x n] is generated by 1/16 times the determinant of 



/y»4 I /-y»4 _|_ <T^4 I ry>4 j " ~" /f»2 I /y»2 ry>2 /y>2 />-»2 I / y* 2 />-»2 ry» 2 /->-» 2 I <-yi 2 /y«2 /y» cy» /y« 

x 00~ l_J '01~ rJ '10~ l ~ J 'll J '00 J '01~ rJ '10 J 'll x 00 J '10~ rJ '01 J 'll J/ 00 J, 11~'~' / '01 X 10 X 00 X 01 J '10 X 11 



Ho 
4< 

Au\ 



2noo«oi 
2ul Q u 01 

2«iou?i 
2uf ttn 



2woi«?i 
2moo«io 
2mqiMii 



2u iU? 

2MqiM10 

2«oo M n 



MoiMioMn 

U 00 ^10^11 
U 00 ^01^11 
MOOMOIMIO 



'IX 



For fixed w^-, this 5 x 5-determinant defines Kummer's quartic as a surface in P 3 , with 
coordinates (xoo : 201 : £io : Xn), written as in [6, Exercise 3, page 204] or [42, page 354]. 

Note that the last four rows in (1.1) represent the Jacobian matrix of the first row. The 
matrix derives from the fact that (uqq : u m : u 10 : Uu) must be a singular point on the 
Kummer surface. The surface and its 16 nodes are invariant under the group (Z/2Z) 4 which 
acts by sign changes and permutations [42, page 353]. The 15 other nodes are easily found: 

(woo : -wio : Uqi : -Wn), («oo : u w : -u 01 : -u n ), (u n : -u 01 : -itio : «oo)- 
A convenient representation of the 166 configuration of nodes is the matrix product 



1.2] 



whose 16 entries are the linear forms whose coefficients are the 16 nodes. As explained in 
Hudson's book [33, §16], the combinatorial structure of the 16 nodes can be read off from this 
4 x 4- matrix, and it leads to various alternate forms of the defining quartic polynomial in [33, 
§19]. In [33, §102] it is shown that the matrix product (1.2) expresses quadratic monomials 
in theta functions with characteristics in terms of the second order theta functions. 
To illustrate how our coordinates can express geometric properties, we note that 
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vanishes if and only if the given abelian surface is a product of two elliptic curves. If this 
happens, then the 16g configuration of all nodes degenerates to a more special matroid, and 
the quartic Kummer surface in P 3 degenerates to a double quadric. The 10 factors in (1.3) 
are the non-zero entries left in the matrix product (1.2) after replacing each Xij by u^. □ 

This article is organized as follows. In Section 2 we review classical material that can be 
mostly found in the books of Coble [11] and Dolgachev-Ortland [18]. We define a Kummer 
threefold in P 7 as the image of a transcendental map whose coordinates are second-order 
theta functions. Each Kummer threefold is the singular locus of the associated Coble quartic 
hypersurface in P 7 , which is a natural genus three analogue of the Kummer surface in P 3 . 

In Section 3 we focus on the moduli space .4.3(2,4) of polarized abelian threefolds with 
suitable level structure. That space is a quotient of the Siegel upper halfspace. It is embedded 
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into P 7 by second-order theta constants [25]. The resulting Satake hypersurface S has degree 
16. Its defining polynomial has 471 terms with integer coefficients, displayed in Proposition 
3.1. We examine how the combinatorics of this polynomial expresses the geometry of .4.3(2, 4), 
notably its hyperelliptic locus, its Torelli boundary, and its Satake boundary. Later in 
Table 3, this is refined to the stratification of S that was found by Glass [26, Theorem 3.1]. 

The theta series in Sections 2 and 3 are defined over the field of complex numbers C, and 
we use floating point approximations for computing them. The resulting algebraic objects, 
however, do not require the complex numbers. This is, of course, well-known to the experts in 
abelian varieties [6]. All the ideals we feature in this paper can be generated by polynomials 
with integer coefficients, and their projective varieties are thus defined over any field. 

In Section 4 we express the Gopel variety Q as the Zariski closure of the image of an 
explicit rational map P 6 -~> P 14 of degree 24. Its 15 coordinates are polynomials of degree 
7 that span an irreducible representation of the Weyl group of type E7. The root system 
and its reflection arrangement play a prominent role, as do configurations of 7 points in P 2 . 
After completion of our work, we learned that our parametrization of the Gopel variety had 
already been studied in [12], under the name Coble linear system. The connection is made 
explicit in Theorem 4.6 where we give an affirmative answer to [12, Question 4.19]. 

In Section 5 we study the defining prime ideal of the Gopel variety Q. We show that it is 
minimally generated by 35 cubics and 35 quartics in 15 variables. This ideal is Gorenstein, 
it has degree 175, and we determine its Hilbert series and its minimal free resolution. 

Section 6 is concerned with a beautiful re-embedding of the Gopel variety Q C P 14 into 
a projective space of dimension 134. The 135 coordinates are the Gopel functions of [18, 
§9.7]. The Weyl group W(E 7 ) acts on these by signed permutations. We construct a toric 
variety T of dimension 35 in P 134 whose intersection with P 14 is Q. Algebraically, the ideal 
of Q is now generated by binomials and linear trinomials. We study the combinatorics of the 
toric ideal of T and its associated convex polytope A. Its 63 distinguished facets and its 135 
vertices are indexed by the lines and the Lagrangians in the finite symplectic space (F 2 ) 6 . 

In Section 7 we return to the embedding of /C3 defined by theta constants and theta 
functions. The universal Coble quartic is an irreducible subvariety of codimension two in 
P 7 x P 7 . It is the complete intersection of the Satake hypersurface S of degree (16, 0) and 
one other hypersurface C of degree (28, 4). The latter is given by an explicit polynomial with 
372060 terms in 16 variables. Together, the two polynomials generate a prime ideal. 

The universal Kummer variety /C3 lives in iSxP 7 C P 7 xP 7 . Other variants of this variety 
can be defined in P 6 x P 7 , via the parametrization of Section 4, or in Qx¥ 7 C P 14 xP 7 , via 
the Gopel embeddings of Sections 5 and 6. In Section 8 we study the bihomogeneous prime 
ideals for these three variants of the universal Kummer variety. We derive lists of minimal 
generators for two of these ideals. Conjectures 8.1 and 8.6 state that these lists are complete. 

Section 9 examines all of our constructions from the perspective of tropical geometry. 

Supplementary materials 

We have created supplementary files so that the reader can reproduce many of the calcula- 
tions that are claimed throughout the text. Most of them are in the format of Macaulay 2 
[27]. These files can be found at http://math.berkeley.edu/~svs/supp/univ_kummer/. 
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A word on the symbolic computations we present: many of the calculations are signifi- 
cantly faster (almost all finishing within a few seconds) if the field of coefficients chosen is 
finite; one that we have used is Z/101. However, we have made no attempt to verify which 
primes give bad reduction. All calculations can also be performed over the rational num- 
bers. In this setting, many calculations take at most a few minutes, and the hardest one in 
Theorem 5.1 took approximately 15 minutes on a modern mid-level performance computer. 
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2 From theta functions to Coble quartics 



In this section, we review Riemann's theta function and its relatives in the case of genus 
3, and we use this to give a definition of Kummer threefolds as subvarieties of P 7 . We also 
discuss background material on the action of the 2-torsion of the associated abelian threefold 
on P 7 and introduce the Coble quartic associated to a non-hyperelliptic Kummer threefold. 

Let r be a symmetric 3x3 matrix with complex entries whose imaginary part is positive 
definite. The set ^3 of all such matrices is a six- dimensional complex manifold, called the 
Siegel upper halfspace. Each matrix t E Sjs determines a lattice A = Z 3 + rZ 3 of rank 6 
in C 3 , and a three-dimensional abelian variety A T = C 3 /A. The Riemann theta function 
corresponding to a matrix r £ ^3 is the function 9 : C 3 — > C defined by the Fourier series 



(r; z) = exp \ [ Txin t Tn + 2irin t z\ . 



(2.1) 



This series converges for all z £ C 3 and r £ #3, and it satisfies the functional equation 



9(r]z + a + rb) = 6(r ; z) ■ exp 



2m(-b t z - -&V&) 



for a, b £ Z 3 



(2.2) 



Deconinck et al. [15] gave a careful convergence analysis and they implemented the numer- 
ical evaluation of 6 in maple. Their work has been extended by Chris Swierczewski who 
implemented the evaluation of abelian functions in Sage [50]. The Riemann theta function 
9 can now be called in Sage as RiemannTheta(tau) (z), where tau is a Riemann matrix and 
z is a complex vector. That code was used extensively in our computations. 

Every pair of binary vectors e, e' £ {0, l} 3 defines a theta function with characteristics 



9[e\e'](r; z) 



ex p 

ngZ 3 



m(n + |)V(n + |) + 27rz(n + ^{z + | 



(2.3) 
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From inspection of this Fourier series, one can see that 

9[e\e'](r;-z) = (-l) eV • 9[e\e'](r; z). (2.4) 

Therefore, of the 2 2 ' 3 = 64 theta functions with characteristics, precisely 2 3_1 (2 3 + 1) = 36 
are even functions of the argument z G C 3 , and the other 2 3_1 (2 3 — 1) = 28 are odd functions 
of z. We shall refer to these as even (or odd) theta functions. 

Finally, for any binary vector a G {0, l} 3 , we consider the second order theta function 



9 2 [a] (r; z) = 9 (2r; 2z + to) ■ exp 



, a l ra , 

m h 2a f z 

2 



(2.5) 

z . 



exp 27ri(n + —J r[n + — ) + 47ri(n + —J 



The second order theta functions are related to the theta functions with characteristics of 
an isogenous abelian threefold by the formula 

& 2 [o-](r;z) = 0[<r|O](2<r;2z). (2.6) 

Further relations between first and second order theta functions are the addition theorem 

e[e\e'}(r; z + w) ■ 6[e\e'}(r; z - w) = ^ (-1) CT ^ 2 [a](r; u>) • 2 [<x + e](r; z) (2.7) 

and its inversion 

8-e 2 [a](T;w)-e 2 [a + e](r;z) = ^ (-1)^' • 9[e\e'](r; z + w) ■ e[e\e'}(r; z - w). (2.8) 

For a fixed matrix r6% the eight second order theta functions define the Kummer map 

k t :C 3 ^F 7 , z ^ (6 2 [000](t;z) : 9 2 [001](r;z) : ■ •• : 2 [111](t; z)). (2.9) 

The identity (2.2) implies that k t factors through a map A T — > P 7 from the abelian threefold 
A T = C 3 /A. This map, which we also denote by k t , has the following geometric description. 
The equation 9{r\ z) — defines the theta divisor on A r . The divisor 20 is ample but not 
very ample. The eight functions (2.5) form a basis for its space of sections. These are known 
as the Schrddinger coordinates, and we denote them by Xoooj^ooi; • • • , x m- The morphism 
A T — > F(H°(A T , 20)) ~ P 7 is given in coordinates by (2.9). The image of the Kummer map 
k t in P 7 is isomorphic to the quotient A T j {z = —z}. We call this the Kummer threefold of r. 

Let A T [2] denote the subgroup of two-torsion points in the abelian threefold A T . This is 
a group of order 64 which we identify with A/2A ~ (F 2 ) 6 . This abelian group acts naturally 
on the set of second order theta functions via 

Q 2 [a](r;z) ^ 6 2 [o-](r; z + 5) where 5 G A T [2]. (2.10) 

This defines an action on P(H°(A r , 20)) ~ P 7 by permuting the coordinates up to sign. 
The action lifts to a linear representation of the Heisenberg group H (see [6, Chapter 6]). 
This mildly non-abelian group is a certain central extension 

1 -> F 2 -> H -> A T [2] -»■ 1. 
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This can be made explicit in terms of the Schrodinger coordinates Xoocb ^ooij • • • > x m we are 
using on P 7 . The Heisenberg group H is generated by the following six operators: 



■Ei,j,k ' ^ x i+l,j,ki x i,j,k ' ^ x i,j+l,ky x i,j,k ' ^ •^i,j,fc+lj ^ 

%i,j,k ' ^ ( 1) ^i,j,ki ^i,j,k ' ^ ( ^i,j,k ' ^ ( 1) ^i,j,k- 

These operators commute up to sign and give the projective representation of A T [2]. 

For any given matrix r G #3, the Kummer threefold has degree 24 in P 7 . When r is 
the period matrix of a smooth non-hyperelliptic curve of genus three, the prime ideal of the 
corresponding Kummer threefold is minimally generated by 8 cubics and 6 quartics in the 
variables x. The Hilbert polynomial of such a Kummer threefold, which is 4n 3 + 4, agrees 
with the Hilbert function for n > 1. This was first derived by Wirtinger in [53, §21]. 

Remark 2.1. In commutative algebra [19], it is customary to also look at higher syzygies. 
Numerical invariants are read off from the Betti table. For a general Kummer threefold, it is 



total 

1 
2 
3 



1 2 3 4 5 

1 14 48 56 24 3 
1 



6 48 56 24 3 



This Betti table shows that the Kummer threefold n T (A T ) is not arithmetically Cohen- 
Macaulay. Failure to be arithmetically Cohen-Macaulay is a general phenomenon for abelian 
varieties of dimension > 2 outside of the case of Kummer surfaces in P 3 (in fact, this is a 
general fact about varieties whose structure sheaf has intermediate co homology via [20, 
Corollary A1.12, Proposition A1.16]). □ 

To compute the above Betti table (in Macaulay 2 [27]), we used the geometric represen- 
tation of the Kummer threefold as the singular locus of a certain quartic hypersurface C T 
in P 7 . That hypersurface was characterized by Arthur Coble in his seminal book [11]. We 
follow Dolgachev and Ortland [18] in our discussion of the Coble quartic C T . 

Proposition 2.2 ([4, Proposition 2.2], [11, §33], [18, §IX.5, Proposition 7]). Let A T be 

the Jacobian of a smooth non-hyperelliptic curve of genus three. There is a unique quartic 
hypersurface C T in P 7 whose singular locus equals the Kummer threefold n T (A T ). The eight 
partial derivatives of the defining polynomial of C r span the space of cubics containing k t (A t ) 
and they generate the prime ideal of n r (A r ) up to saturation. 

The quartic C T is invariant under the action of the Heisenberg group H. The space of 
iJ-invariant quartics is 15-dimensional, and the defining polynomial F T of the Coble quartic 
hypersurface can be written as a linear combination of a basis for this space of if-invariants. 
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Using the Schrodinger coordinates x on P 7 , 



we write 



<Y> , ( ' rf I iy>^ I fy»4 I ^y«4 I nr' I T* I nr I rf ^ 

' l x ooo ^ -^ooi ^ •''Oio t x 011 -|- x 100 -|- x 101 -|- x 110 -t- X ln ) 




I q _ „ _ . /" / y* 2 fy»2 I , y» 2 fy»2 I />-»2 ,y»2 I ,-y-» 2 ,-v-»2 \ 
-|- iiii l J -000 x lll "r x 001 x 110 x 010 x 101 x oii J aooJ 

+ ^001 ' (^000^010^100^110 + ^OOl^Oll^lOl^lll) 

+ ^oio ' (^000^001^100^101 +^010^011^110^111) 
+ ton ' (^000^011^100^111 +^001^010^101^110) 
+ ^100 ' (^000^001^010^011 +^100^101^110^111) 
+ tioi • (^OOO^OlO^lOl^lll + ^ooi^on^ioo^no) 
+ tno • (^OOO^OOl^llO^lll + ^OlO^Oll^lOO^lOl) 

+ tm • (xooo^on^ioi^no + ^001^010^100^111)- 



(2.12) 



This representation appears in [18, §IX.5, Proposition 8] and after Theorem 3.2 in [5, Section 
3]. The 15 coefficients r, s,,t, are parameters. This notation is used throughout this paper. 

Remark 2.3. The monomials in x that appear in the equation (2.12) of the Coble quartic 
can be understood combinatorially via the affine geometry of the eight-point vector space 
(F 2 ) 3 . Namely, the four terms multiplied by are the four affine lines parallel to the vector 
(i, j, k), and the two terms multiplied by are the two affine planes perpendicular to the 
vector (z, j, k), with respect to the usual dot product in (F 2 ) 3 . □ 

In this paper, each non-hyperelliptic Kummer threefold will be represented as the variety 
in P 7 cut out by the eight partial derivatives dF T /dxijk- The adjective "universal" in our 
title means that we are working over the six-dimensional base of all Coble polynomials F T . 
This family is obtained by letting r run over the Siegel upper halfspace #3. In particular, 
r, s.,t. depend analytically on r. We shall review the relevant moduli space in Section 3 

The Coble quartic is closely related to the Kummer surface in Example 1.1. Indeed, the 
Kummer quartic is the expansion of the determinant (1.1) along the first row: 



The monomials in x 0) ^oi) ^ioj ^11 seen here correspond to the affine subspaces of (F 2 ) 2 . The 
coefficients are polynomials of degree 12 in the theta constants u^. They are obtained as 
the 4 x 4-minors of the last four rows in (1.1). These minors satisfy the cubic equation 



This cubic defines a hypersurface in P 4 which is known as Segre 's primal cubic. In Section 5 
we shall derive the analogous relations for the fifteen coefficients of (2.12). 



r • (xoo+^oi+^Io+^ii) + * ' (^00^01^10^11) + 
s 01 ' (^oo^oi+^io^ii) + S 10 ' (^OO^lO+^Ol 2 -!!) + s n ' (^oo^n+^oi^io)- 



(2.13) 



16r 3 + rt 2 + 4(s iSioSii - rs 2 



01 



rs 2 10 -rs 2 n ) = 0. 



(2.14) 
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Remark 2.4. Both the Kummer surface in P 3 and the Coble quartic in P 7 are self-dual 
hypersurfaces [45]. The role of the 16 singular points on the Kummer surface is now played 
by the 64 singular points on the Kummer threefold. These are the images of the 2-torsion 
points of A T under the Kummer map k t . In analogy to the 16 entries of the matrix (1.2), we 
consider the 64 linear forms in (3.10) below. Their coefficients give the 64 special singular 
points of k t (A t ) in P 7 . These points lie on 64 special hyperplanes. Self-duality gives a 64 2 s 
configuration consisting of points and hyperplanes in P 7 . If A T is the Jacobian of a genus 
three hyperelliptic curve, then the Coble quartic F T becomes the square of a quadric, as we 
shall see in (8.9), and the Kummer threefold fails to be projectively normal [35, §2.9.3]. □ 

3 The Satake hypersurface 

In order to understand the family of all Kummer threefolds, we now vary the matrix r 
throughout the Siegel upper halfspace ^3. The modular group Sp 6 (Z) consists of block 

matrices 7 = where a, b, c, d are 3x3 matrices with integer entries such that 7 Jrf = 

J, where J = 0y ' following [6, 18, 42], the modular group Sp 6 (Z) acts on fj 3 by 

7 or = (rc + d)~ l (to + b) . 
The quotient is the moduli space of principally polarized abelian threefolds: 

^3 = iVSp 6 (Z). 

We will also consider certain level covers of .A3, which can be constructed by taking 
quotients of fj 3 by appropriate normal congruence subgroups of Sp 6 (Z). The subgroup 

T 3 (2) = ( 7 e Sp 6 (Z) : 7 = ( U n 3 T ? ) (mod 2)1 (3.1) 



x M 3, 

has index |Sp 6 (F 2 )| = 1451520; see (6.2). We also define 

T 3 (2,4) = {7=^ Jjer 3 (2) : diag(a T 6) = diag(c T c/) = mod 4 

This group has index 64 in T 3 (2), and the quotient group r 3 (2)/r 3 (2, 4) is isomorphic to 
(F 2 ) 6 . These subgroups determine moduli spaces 

^ 3 (2) = iVr 3 (2) and A 3 (2,4) =i3 3 /r 3 (2,4). 

We record that the induced quotient map of moduli spaces is a 64-to-l cover: 

^3(2,4) ^>A 3 (2). (3.2) 

The Torelli map gives an embedding of the moduli space A^ 3 of smooth genus three curves 
into A3. We can thus regard .M3 as a subset of A3. (For the experts: this would not be 
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correct if we considered these spaces as stacks, but we will not do this here.) The inverse 
images of A4 3 in -4.3(2) and in -4 3 (2,4) are denoted by -M 3 (2) and A4 3 (2,4) respectively. 

The moduli spaces above are six- dimensional quasi-projective varieties, and we are inter- 
ested in the homogeneous prime ideals of certain embeddings. In this section we focus on 
.4.3(2,4). Our point of departure is the theta constant map : ^3 — >■ P 7 , which is defined by 



are modular forms of weight 1/2 with respect to the subgroup 1^(2,4) of Sp 6 (Z). Hence 
the theta constant map d can be regarded as a morphism from the level cover .43(2,4) = 
fj 3 /r 3 (2,4) into the projective space P 7 . This map is an embedding by [26, Theorem 4.1]. 

The closure S = $(.4.3(2,4)) of the image of the theta map in P 7 is a six- dimensional 
hypersurface. It is isomorphic to the Satake compactification „4 3 (2,4) of the moduli space 
.4 3 (2, 4). This is a degree 16 hypersurface which we call the Satake hypersurface. Its defining 
polynomial, which we also denote by S, was found by Runge [46, §6] and we now describe it. 
Write G for the permutation group of order 1344 that is generated by the four transpositions 

(%oi«oio)(«ioi«iio), (uoiou 100 )(u 01 iu 101 ), {u 00 ou Q01 )(u ow u 011 ), and («ioo«ioi)(«iio«ni)- 

As an abstract group, G = SL 3 (F 2 ) x (F 2 ) d . In our context, it is precisely the subgroup 
of Sp 6 (F 2 ) that acts on the Schrodinger variables u by coordinate permutations. Let [u*] m 
denote the sum over the G-orbit of a monomial u*. The index m is the size of this orbit. 

Proposition 3.1 (Runge). The Satake hypersurface S = $(„4 3 (2, 4)) is an irreducible hy- 
persurface of degree 16 in P 7 . Its defining polynomial is the following sum of 471 monomials: 

[Mq 00 Mq 01 Mq 10 M^ 00 ] 56 — 2[Mq 00 M 001 M 010 M 1 ooMo11 m 101 m 110 m 11i]8 + 2[Wooo M 001 M 010 M 100 M 01l' U 10l]84 
+ K)0 u 001 u 010 u 100 u lll]56 - [MooO M 001 M 010 M ?00 M llo]224 + 4[M^ 00 M^ 01 M010Ml00M01lMl0lMll0Mlll]28 

- 16 H 00 ul 01 u 3 ow u W ou 3 ou u wl u lw u lu ]u + 72ul 00 ul 01 ul 10 ul 00 ul 11 u 2 101 ul 10 uj 11 . 

A more compact representation of the Satake hypersurface S, in terms of second order 
theta constants, appears in [25, Example 1.4]. Our 471 monomials can be derived from these. 

We next discuss a beautiful direct relationship between the polynomials in Example 1.1 
and Proposition 3.1. We learned this from Sam Grushevsky and Riccardo Salvati Manni. 
Both polynomials have total degree 16 in eight unknowns, and the former appears as an 
initial form in the latter. Let q denote a deformation parameter and set 

uooo = u 00 + 0(g 4 ), u 001 = u m + 0(g 4 ), u 010 = u w + 0(g 4 ), u n = u n + 0(g 4 ), ^ g . 
U100 = 2qx 00 +O{q 9 ), u wl = 2qx 01 +O(q 9 ), u lw = 2qx w +0{q 9 ), u in = 2qx u +0{q 9 ). 

Under this substitution, the Satake polynomial in Proposition 3.1 takes the form 



r M. (O 2 [000](r;0) : O 2 [001](r;0) : ••• : e 2 [lll](r;0)). 



(3.3) 



This map is not injective, since the second order theta constants 



O 2 [a](r;0) = ^ exp 2<Ki(n + |)V(n + |) 



(3.4) 



s 



det[.]-g 4 + 0(g 8 ), 
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where [•] is the 5 x 5-matrix of Example 1 . 1 that defines the universal Kummer surface /C2 • 
This identity can be derived from the Fourier- J acobi expansion of second-order theta 
constants. Namely, in that expansion we write the 3 x 3-matrix r G S) 3 in the form 



r = ( ^ T i J 1 where r' G #2, s G C, and z G C 2 . 
If we set q = e 7 ™ 5 / 2 , then each of the eight genus three second order theta constants 




u a = 6 2 (r;0)[a] 



has a Taylor series expansion in q. The leading coefficient in these Taylor series is either a 
theta constant or a theta function of genus two. These expansions are given in (3.5). 

The Satake hypersurface S = $(.4.3(2, 4)) in P 7 contains several loci of geometric interest: 

• The hyperelliptic locus has codimension one in S, and hence codimension two in P 7 . 

• The Torelli boundary S \ $(.M3(2,4)) has codimension two in S. 

• The Satake boundary S \ $(.4.3(2, 4)) has codimension three in S. 

For each of these loci, we shall describe its irreducible components and defining polynomials. 
We begin with the Satake boundary. By [24, Lemma 3.5], it consists of 126 three-dimensional 
subspaces P 3 in P 7 : each of the 63 non-zero half-periods e G A/2A induces a linear involution 
on P 7 via the action on second order theta constants given in (2.10), and the fixed point set of 
this involution on P 7 is the union of two P 3 s. For instance, for e = (1/2, 0, 0), the involution 
(2.10) fixes the coordinates u oo, u ooi, u ow, uou, it switches the sign on Uwo, wion uuo, uiu, 
and the two P 3 s are obtained by setting either of these two groups of four variables to zero. 

The hyperelliptic locus in S is the closure of the set of all points i?(r) where A T is the 
Jacobian of a smooth hyperelliptic curve of genus 3. It is known (see e.g. [26]) that a genus 
3 curve is hyperelliptic if and only if one of its 36 first order theta constants #[e|e'](r;0) 
vanishes. We write these 36 divisors in our eight coordinates mqoo, • • • , um using the formula 



which is obtained by setting w = z = in the addition theorem (2.7). Hence the hyperelliptic 
locus in S can be defined set-theoretically by the equation of degree 72 obtained by taking 
the product of the quadrics (3.6) where (e, e') runs over the 36 even theta characteristics. 

We will see in Section 7 that the product Y[ $[ e l e ']( r i 0) of the 36 first order theta constants 
is in fact a polynomial in u, so the scheme defined by the degree 72 polynomial is not reduced. 
Modulo each of the quadrics (3.6), the Satake polynomial becomes the square of an octic. 
In the supplementary materials, we give this octic for one of the 36 components, as well 
as a verification that the subscheme of this component defined by the octic is reduced. So 
each component of the hyperelliptic locus is a complete intersection of degree 16. Hence as 
a subvariety of P 7 , the hyperelliptic locus has codimension two, degree 576 = 16 ■ 36, and 
36 irreducible components. Glass analyzed in [26, Theorem 3.1] how many of these 36 theta 
constants can simultaneously vanish. We shall refine his results in our Table 3. 

The Torelli boundary is the closure of the image under $ of the set of all polarized abelian 
threefolds that decompose as the product, as a polarized variety, of an elliptic curve and an 
abelian surface. This means that r can be transformed into block form under Sp 6 (Z). 




(3.6) 
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Proposition 3.2. The Torelli boundary coincides with the singular locus of the Satake hy- 
persurface S. It is the union of 336 irreducible four- dimensional subvarieties ofF 7 , each of 
which is defined by the 2 x 2-minors of a 2 x 4-matrix of linear forms in Mooch M ooi) • • • ? u m- 

Proof. By [25, Lemma 3.2], the Satake hypersurface is non-singular at each point that repre- 
sents an indecomposable polarized abelian threefold. It follows from [26, Theorem 3.1] that 
a polarized abelian threefold decomposes as a product if and only if at least two of its first 
order theta constants vanish. Moreover, the vanishing of two first theta constants implies 
the vanishing of at least six of the quadrics in (3.6). The relevant 6-tuples I of even theta 
characteristics m = [e|e'] have the property that, for any three m 1 ,m 2 ,m 3 in J, the sum 
mi + m 2 + 777,3 is an odd characteristic. There are 336 such 6-tuples /. They correspond 
to azygetic triads of Steiner complexes, and hence to non-isotropic planes in the symplectic 
vector space (F 2 ) 6 , and hence also to root subsystems A 2 in E 7 ; see [39, Prop. 1(2)]. 

Now, if $(t) is in the Torelli boundary then the 3 x 3- matrix r is in the Sp 6 (Z)-orbit 
of a matrix tq that decomposes into two blocks given by matrices in fji and .fj 2 . Since 
theta constants behave multiplicatively under this decomposition, the image of the locus of 
To admitting such a decomposition is a Segre variety P 1 x P 3 in P 7 . That Segre variety is 
defined by the 2 x 2-minors of the 2 x 4-matrix 

/ TfOOO 7^001 U 010 u 01l 
V u 100 UlOl "HO "111 

It is easy to check in Macaulay 2 that the ideal of 2x2 minors of (3.7) contains the partial 
derivatives of the Satake polynomial S. In particular, this component consists of singular 
points. Since the other 335 components P 1 x P 3 are obtained by applying the action of the 
modular group, we see that all components are in the singular locus. Putting everything 
together, we conclude that the singular locus of S coincides with the Torelli boundary. □ 

The 336 irreducible components of the Torelli boundary in S are a direct generalization 
of the 10 irreducible factors in (1.3). In both cases, the components are defined by certain 
quadrics in the u-coordinates that can be written as 2 x 2-determinants of linear forms. 
In Section 7 we shall return to the equations of the Torelli boundary when we study the 
universal Coble quartic in P 7 x P 7 . Let us now define our universal objects in precise terms. 

We combine the Kummer map (2.9) and the theta constant map (3.3) by setting 

k: S) 3 x C 3 ^P 7 x P 7 , (t,z) ^ ($(t),k t (z)). (3.8) 

This is the universal Kummer map in genus g = 3. The closure of its image in P 7 x P 7 is 
the universal Kummer threefold. This irreducible variety of dimension nine is denoted 

£ 3 = £ 3 (2, 4) := k(£ 3 x C 3 ). 

We fix the coordinates u = (u 000 : u ooi Um) on the first copy of P 7 , and we fix 

the coordinates x = (x oo : ^001 : ■ ■ • : ^m) on the second P 7 . Thus the represent 
the theta constants which can be regarded as coordinates for the moduli space ^(2,4) of 
polarized abelian threefolds with level structure, while the represent the second order 
theta functions which are coordinates for the individual abelian threefolds A T themselves. 
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The following two formulas relate these coordinates on our P 7 x P 7 to the theta functions 
with characteristics. First, by specializing w = z in (2.7) we obtain the formula 

9[e\e'](r;2z)-9[e\e'](r;0) = £ (-1)™' • x a ■ x a+e . (3.9) 

Second, the specialization w = in (2.7) yields 

9[e\e'}(r;z) 2 = ^ (-1)^ • u a ■ x a+e . (3.10) 

The identity (3.6) arises from either of these by setting z = 0. Our formulas admit Z[l/8]- 
linear inversions, derived from (2.8), that express the various quadratic monomials on P 7 x P 7 
in terms of thetas with characteristics. Equation (3.9) expresses first order theta functions 
with doubled argument 2z as quadratic polynomials in the second order theta functions. 

It is sometimes advantageous to embed the Kummer threefold not in P 7 , but in the larger 
space P 35 whose coordinates are indexed by even pairs (e, e'). Likewise, it makes sense to re- 
embed the universal Kummer variety /C 3 (2, 4) from P 7 x P 7 into P 35 x P 35 . Using the addition 
theorem for theta functions (2.7), this can be accomplished in two different ways. First, we 
can use the formulas (3.10) and (3.6) to map £3(2,4) into P 35 x P 35 with coordinates 

(9[e\e'}(r;0) 2 ,9[e\e'}(T;zY). (3.11) 
Second, we can use the formulas (3.9) and (3.6) to map £3(2, 4) into p 35 xP 35 with coordinates 

(9[e\e'}(r;0) 2 ,9[e\e'](r;2z)). (3.12) 
In both cases, the re-embedding is given by a certain Veronese map P 7 x P 7 — > P 35 x P 35 . 
Remark 3.3. Our motivating problem was to determine the bi-homogeneous prime ideal 

X 3 C Q[u,x] 

of the universal Kummer threefold £3(2,4). One of the minimal generators of X 3 is the 
Satake polynomial of degree (16,0). To find others, one might try the following approach. 
For any two non- negative integers r and s, consider the space Q[u, x] r s of polynomials that 
are bi-homogeneous of degree (r, s). This space has dimension ( ^ r ) ( 7 ^ s ) ■ We seek to identify 
the subspace (X 3 )( riS ) of polynomials that lie in our ideal X 3 . That subspace can be computed 
using (numerical) linear algebra. The basic idea is simple: using Swierczewski's code for the 
Riemann theta function 9, we implemented pieces of Sage code for the second order theta 
functions (2.5), for the Kummer map (2.9), for the theta constant map (3.3), and for the 
universal Kummer map (3.8). For each point (r,z) E S) 3 x C 3 , we can thus compute one 
linear constraint on polynomials in (X 3 )( r s ), as these vanish at k(t, z). By plugging in enough 
points, we get linear equations in ( 7 ^ r ) ( 7 ^ s ) unknowns whose solution space equals (X 3 )( r s ). 

In practice, however, this approach does not work at all. The primary reason is that 
the size of our linear systems is too large, even with the use of state-of-the-art software for 
numerical linear algebra. The key to any success would be the identification of subspaces in 
which equations for X 3 can lie. Methods from representation theory are essential here. 
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For instance, suppose we are told that the Satake hypersurface has degree 16, and we 
are asked to find its polynomial S. Solving the naive system of linear equations in ( t, 16 ) = 
245157 unknowns is an impossible task. However, all but 471 of the unknown coefficients are 
zero. Solving the restricted system for the 471 coefficients is easy. We do this by computing 
a singular value decomposition of a matrix whose rows are the 471 monomials evaluated at 
$(Tj) for sufficiently many 3 x 3-matrices Tj. The same works with the 1168 monomials of 
degree (16,4) in Lemma 8.2. From these we can compute a polynomial in (I3)(i6,4)- Q 



4 Parametrization of the Gopel variety 

Consider the family of homogeneous quartic polynomials in eight unknowns given by 
the formula in (2.12). This family constitutes a projective space P 14 whose coordinates are 

(r : Sooi : Soio : ' ' ' : s m : ^ooi : ^oio : ' • " : inij- 

The Gopel variety is the six-dimensional irreducible subvariety Q of P 14 that is obtained as 
the closure of the set of all Coble quartics F T where r runs over non-hyperelliptic r £ S)%. 

Coble [11, §49] describes 63 cubic polynomials that cut out Q set-theoretically; see also 
[18, §IX.7, Proposition 9]. In Section 5 we review Coble's construction, and we determine 
the prime ideal and its graded Betti table. The Gopel variety Q is a compactification of the 
non-hyperelliptic locus inside the moduli space ^3(2). See [37] for an analytic approach. 

We now present an utterly explicit polynomial parametrization 7 of the Gopel variety Q. 
Let P 6 be the projective space with coordinates [c\ : C2 : C3 : C4 : C5 : Cq : oj). Our map 

is defined by the following 15 homogeneous polynomials of degree 7 in ci, c 2 , . . . , c 7 : 



r = Ac^CsC^CqCj 



•S001 


= cic 2 c 7 (c 4 


2c 3 C 4 -\- C 4 


- 2c 2 c 2 


- 2c 2 c 2 + c 4 


- 2r 2 r 2 


- 2r 2 r 2 


- 2c 2 c 2 + 4) 


s 010 


t 4 

— — ClC3 c 5\ c 2 


— 2c 2 C 4 -\- C 4 


- 2c 2 c 2 


_ 9r 2 r 2 + r 4 




2C^Cy 


2cgC-^ "t~ 


s 011 


1 4 

— CiC4Cq{C 2 


— 2c 2 c 3 + c 3 


- 2r 2 r 2 


- 2c 2 c 2 + 4 




2C^Cy 




SlOO 


= -C 2 C3C 6 (Ci 


- 2c-^c 4 -\- c 4 


- 2c 2 c 2 


~ ?44 + 4 




2C4C7 


2CgC^ - r~ CyJ 


SlOl 


( 4 

— C2 C A C 5\ C \ 


— 2c 1 c 3 + c 3 


- 2c 2 c 2 


_ 9r 2 r 2 + r 4 




2C^Cy 


2cgC^ ~\~ Cij^j 


SllO 


= -C2,C A C 7 (c{ 


— 2c 1 c 2 + c 2 


- 2c 2 c 2 


— 2c 2 c 5 + c 5 


- 2c 2 c 2 




- ^14 + 4) 


s lll 


( 4 

— C5 c 6 c 7\ c l 


— Ic^c^ + c 2 


- 2c 2 c 2 


— 2c 2 c 3 + c 3 


2c^c 4 


— 2c 2 c 4 


— Qc^c^ -\- c 4 ) 


^100 


i 4 2 

= ci(-c 2 c 3 - 


f- c 2 c 3 + c 2 c 4 


4 2 
c 3 c 4 


2 4,24 

- c 2 c 4 -+- c 3 c 4 


4 2 
C 2 C 5 


■ 2r 2 r 2 r 2 


+ 2c 2 c 2 c 2 + c. 



1 24_ 24, 42_ 42_|_9222_9222_|_9222_9222_ 42 

' C 2 C 5 C 4 C 5 C 2 C 6 C 3 C 6 ZC 2 C 4 C 6 /C 3 C 4 C 6 /C 3 C 5 C 6 /C 4 C 5 C 6 C 5 C 6 

2 4 , 2 4 , 2 4 , 9 2 2 2 , 4 2 9 2 2 2 4 2 , 9 2 2 2 

C 2 C 6 ' C 3 C 6 ' C 5 C 6 ' ZC 2 C 3 C 7 1 C 3 C 7 ZC 2 C 4 C 7 C 4 C 7 -f- ZC 2 C 5 C 7 

_ 9 2 2 2 , 42_9222|_9222_ 42_ 24, 24_ 24, 2 4\ 

, I 42,24,42 42 24,24,42 4 2 , r> 2 2 2 r> 2 2 2 

t-010 — c 2\ — c l c 3 ~r C^C 3 -+- C 4 C 4 — C 3 C 4 — C^C 4 + C 3 C 4 + C^Cg — c 3 c 5 + zc 4 c 4 c 5 — zc 3 c 4 c 5 

24, 24_ 42_9222_ | _9222_ | _ 42_,_9222_9222 
C 1 C 5 1" C 3 C 5 C 1 C 6 ZC 1 C 3 C 6 1" /C 3 C 4 C 6 + C 4 C 6 + /C 3 C 5 C 6 /C 4 C 5 C 6 
, 42, 24_ 24_ 2419222, 42_9222_ 42_9222 
' C 5 C 6 ' C 1 C 6 C 4 C 6 C 5 C 6 + ZC 1 C 3 C 7 ' C 3 C 7 ^C^C^ C 4 C 7 ZC 1 C 5 C 7 

_i_9222_ 42_|_9222_9222_,_ 42_ 2 4, 2 4, 24_ 24\ 
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i / 42,24,42 42 24,24,0222,42 r>222 

^OOl — — CjC 2 + C^C 2 "T C^C 4 — C 2 C 4 — C-^C^ -+- C 2 C 4 ■+- ZC-^^Cg ■+- C 2 Cg — ZC^C^C^ 

42_ 24, 24_ 42_9222|_9222, 42,r)222_9222 

42, 24_ 24, 24, 42_ 42_|_9222_9222_9222 

,9222, 42|_9222_9222_ 42_ 24, 24_ 24, 2 4\ 
"I - ZC 4 Cr^C 7 "t - Cr^Cj \~ ZC^CqC^ ZC 4 CgC 7 CgC 7 Ci C 7 "t - C 2 C 7 Cr^Cy "t - CgC 7 J 

. l'_ 42 i 2 4 , 4 2 _ 42 _ 2 4 , 24 _ 42_9222,9222_,_ 42 

till — c 4{ c l c 2 1* CiC 2 -+- C1C3 C 2 C 3 C1C3 -+- C 2 C 3 C^Cg ZCiC 2 C 5 -+- ZC 2 C 3 C 5 -|- C 3 C 5 

, 24_ 24_|_9222, 42_9222_ 42_|_9222_9222, 42 

-+- C 1 C 5 C 3 C 5 -+- ZC 1 C 2 Cg -+- C 2 C 6 /^CgCg c 3 c 6 -+- zc 1 c 5 c 6 /c 2 c 5 c 6 -+- c 5 c 6 

24, 24_ 24, 42_ 42|9222_9222_,_9222_9222 
C 2 C 6 ~r C 3 C 6 C 5 C 6 ' c l c 7 C 2 C 7 1 ZC^C^Cj ZC 2 C 3 C 7 + ZC 2 C 5 C 7 ZC 3 C 5 C 7 

42_9222,9222_,_ 42_ 24, 24, 24_ 2 4\ 

« /_ 42, 24_i9222, 4 2 _ 2 4 _ 4 2 _ 9 2 2 2 _, 9 2 2 2 _ 9 2 2 2 

Oil — ^5 v ^1^2 ~~ ^1^2 ~~ 1^2^3 ~~ ^2^3 ^2^3 ^1^4 ^1 ^2 ^4 ~~ ^C-^C^C^ ^C2C^C^ 
42, 24, 24, 42_ 42_9222, 42,9222_ 42 

C 3 C 4 + CiC 4 + C 3 C 4 -+- CiCg C 2 Cg ZCiC 3 Cg -+- c 3 c 6 -+- zc 2 c 4 c 6 c 4 c 6 

24, 24_ 24, 24, 42_ 42_9222, 42,9222 
C-^ Cg ~r~ ^2^6 ^3^6 ~~ ^4^6 ~~ ^1^7 ^2^7 ^C^C^C^ — t - C3C7 — r~ ZiC^C^C'j 

4 2 I 9 2 2 2 _n 2 2 2 , n 2 2 2 _n 2 2 2 _ 2 4, 2 4 _ 2 4 _i_ 2 4> \ 
. / 42,24 42 0222,24,0222,42,0222 

*ioi = c 6 (-CiC 2 + CiC 2 - CiC 3 - 2c 1 c 2 c 3 + C1C3 + Zc^c^ + c 2 c 4 + 2ciC 3 c 4 

_9222, 42_ 24_ 24, 42_ 4 2 , n 2 2 2 _ 42_9222, 42 

— zc 2 c 3 c 4 -+- c 3 c 4 c 2 c 4 c 3 c 4 -+- c^^Cg c 2 c 5 -+- zc 2 c 3 c 5 c 3 c 5 ZC^C^C^ -+- c 4 c 5 

24, 24, 24_ 24, 42_ 42,9222_ 42_9222, 42 
C-^C^ ~\~ C 2 Ce^ "t - C 3 Ctj ^4^5 ~~ f~ ^1^7 ^2^7 2^3^7 ^3^7 ^C-^C 4 C 7 — p C 4 C 7 

i9222_9222_9222_,_9222_ 24, 24, 24_ 2 4\ 

~\~ ZC-^C^C*j ZC^C^Cj ZC^Cr^Cy ~T~ ZC^C^Cj C-^C'j \~ C 2 C 7 \~ C^Cf C^Cy J 

i i 42,9222 42,24,24 42,9222 42 9222 

6iiq — ^1^3 ' ^1^2^3 ^2^3 ' ^1^3 ' ^2^3 ^1^4 1^2^4 ^2^4 1^3^4 

_9222, 24, 24, 42_9222, 42_ 42,9222_ 42_ 24 
ZC 2 C 3 C 4 -+- C-^C^ -+- C 2 C 4 -+- CiC 5 ZC1C2C5 -+- C 2 C 5 C 3 C 5 -+- ZC 3 C 4 C 5 C 4 C 5 C1C5 

24, 24, 24, 42_9222, 42_ 42,9222_ 42,9222 
C 2 C 5 1 C 3 C 5 + C 4 C 5 + C 1 C 6 /C 1 C 2 C 6 ' C 2 C 6 C 3 C 6 ' /C 3 C 4 C 6 C 4 C 6 ' /C 1 C 5 C 6 
i_9222_9222_9222_ 24_ 24, 24, 2 4\ 
T ZC 2 C 5 C 6 /C 3 C 5 C 6 /C 4 C 5 C 6 C 1 C 6 C 2 C 6 ' C 3 C 6 + C 4 C 6/ 

Theorem 4.1. T/ie Gopel variety Q = As(2) is the closure of the image of the rational map 
7: P 6 P 14 given by the polynomials above. This map defines a 24-to-l cover of Q. 

The rest of this section is devoted to a conceptual derivation and geometric explanation 
of the map 7. This will then furnish the proof of Theorem 4.1. In Section 6 we shall see that 
7 is equivalent to the embedding via Gopel functions described in [18, §IX.8, Theorem 5]. 
Our construction here is based on the representation-theoretic approach developed in [29]. 

Consider the Heisenberg group H that is generated by the six operators in (2.11). Let 
H' be the subgroup of GL 8 (C) generated by H and the scalar matrices C*. Let N(H') be 
its normalizer. Then N(H')/H' = Sp 6 (F2) by [6, Exercise 6.14]. Since every element of 
H' acts on the 15 parenthesized polynomials in (2.12) by scalar multiplication, the space 
spanned by these is preserved by the action of N(H'). This action factors through H, so 
we get a 15-dimensional representation of N(H')/H. This contains Sp 6 (F 2 ) as a subgroup, 
and this 15-dimensional representation of Sp 6 (F 2 ) is irreducible. There is a unique such 
irreducible representation, so we call it U 15 . For our computations we shall use the larger 
group ^(£7) = Sp 6 (F 2 ) x {±1}. This is the Weyl group of the root system of type E 7 . We 
use the same symbol U 15 to denote the corresponding irreducible representation of W^(E 7 ). 

Let t) = C 7 be the reflection representation of W(E 7 ). We will work in two different 
coordinate systems of h. The first one we call c\, . . . , c 7 (with the standard quadratic form 
Si c ?)' which has the property that each q = is a reflection hyperplane. Explicitly, the 
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following two matrices are a pair of generators for W"(E 7 ) in terms of this basis: 
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(4.2) 



One can verify with GAP [22] that they generate a group of the correct size, and it contains 
seven reflections that satisfy the Coxeter relations of type E7. See our supplementary files. 

All representations of the Weyl group jy(E 7 ) are isomorphic to their duals, as is true for 
any real reflection group. In what follows we shall distinguish between representations and 
their duals to avoid confusion. The coordinate ring of f) is Sym(f)*) = C[ci, c 2 , c 3 , c 4 , c 5 , c 6 , c 7 ]. 
There is a unique copy of the 15-dimensional irreducible representation U 15 in Sym 7 ([)*). We 
call this subspace t/ c 15 . It can be realized explicitly by taking the ty(E 7 )-submodule spanned 
by the monomial c^c^c^CqC^ . This is an example of a Macdonald representation [38]. 

Since U 15 = (U 15 )*, there exists a unique basis in £7 15 which is dual to the basis of the 
15 Heisenberg-invariant quartics of (2.12) in U lb . We computed that dual basis of £7 C 15 . It 
was found to consist of the above 15 polynomials r, sooi, • • • , tm of degree 7 in ci, C2, . . . , c 7 . 
These elements of Sym 7 (f)*) are the coordinates of our rational map 7: P 6 ---> P 14 . 

By multiplying the dual bases and summing, we get a jy(E 7 )-invariant in U 15 £g> U^ 5 . 
This is the defining polynomial F T of the Coble quartic, with coefficients as above. This 
concludes our derivation of the map 7. Here is an alternative way of getting our formulas. 

Remark 4.2. Let A be an 8-dimensional complex vector space with basis 

a-i = ^ooo; 02 = X100, a>3 — ^oio ; a i — ^noj «5 = ^ooi ; o-e = ^ioi> a-j = xon, as = xm. 

The Heisenberg group H acts on /\ 4 A and the space of invariants has the following basis: 

h\ = a± A a2 A 03 A 04 + as A a§ A a 7 A a%, 

h>2 = «i A <22 A 05 A a$ + as A 04 A a 7 A as, 

/13 = ai A 03 A 05 A o 7 + a2 A 04 A a§ A a§, 

/14 = a\ A 04 A a§ A a 7 + ai A 03 A 05 A as, 

fi5 = ai A a3 A ag A as + a2 A 04 A «5 A a 7 , 

/i6 = ai A 04 A 05 A as + 02 A 03 A a$ A a 7 , 

/i 7 = ai A 02 A a 7 A a§ + 03 A 04 A 05 A a§. 

The techniques in [29, §6.2, §6.5] show how to start with any vector in /\ 4 A and produce 
the equations for a Kummer variety in the projective space of lines in A* along with the 
equation for its Coble quartic. If we apply these techniques to c±h\ H — • + c 7 /i 7 , then we get 
the above formula for F T . The derivation of Conjecture 8.1 will discuss this in detail. □ 

We now come to the proof that our parametrization of the Gopel variety is correct. 
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Proof of Theorem 4-1- The Macdonald representation f/ 15 in Sym 7 (f)*) is unique. It defines 
a unique W / (E 7 )-equivariant rational map 

7: P 6 — > P 14 . 

Hence, up to choosing coordinates on the two projective spaces, the map 7 agrees with the 
one described in [18, §IX.7, Remark 7]. The closure of its image must be the Gopel variety Q. 

To see this more explicitly, we now come to our second coordinate system, denoted 
d = (di,d 2 , d 3 , g?4, d 5 , d 6 , d 7 ), for the reflection representation f) of W(E 7 ). It is defined by 



Cl 


= d 2 


+ d 4 


+ d 5 


c 2 


= di 


+ d 4 


+ d 7 


c 3 


= d 2 


+ d 3 


+ d 7 


C4 


= di 


+ d 3 


+ d 5 


c 5 


= d 1 


+ d 2 


+ d 6 


c 6 


= d 5 


+ d 6 


+ d 7 


c 7 


= d 3 


+ d 4 


+ d 6 



(4.3) 



The labels on the right hand side define a Fano configuration like (6.3). In these coordinates, 
the 63 reflection hyperplanes of W (E 7 ) are given by the 21 linear forms dj — dj, the 35 linear 
forms dj + dj + d&, and the 7 linear forms d, + dj + d^ + + d m + d re . Here the indices are 
distinct. As a warning to the reader, we note that the d variables do not represent orthogonal 
coordinates for E 7 since the quadratic form £\ cf is not mapped to df . 

The Weyl group W(E 7 ) acts by Cremona transformations on the space of configurations 
of seven points in P 2 . Following Dolgachev and Ortland [18, §IX.7, Remark 7], there exists a 
canonical iy(E7)-equivariant birational map from P(h) = P 6 to that space. In coordinates, 
this canonical map takes (di : d 2 : • • • : d 7 ) to the column configuration of the matrix 

/l 1 1 1 1 1 1\ 

D — I di d 2 d 3 d 4 d 5 d 6 d 7 . (4.4) 
\df d^ d 3 3 d\ d\ d% d 3 7 J 

The rational map 7 is the composition of that birational morphism and the Gopel map of 
[18, §IX.7, Thm. 5], written explicitly in Section 6, that takes 7-tuples in P 2 to points in Q. 

Our claim about 7 being a 24-to-l cover follows from the fact that a general net of cubics 
in P 2 contains 24 cuspidal cubics [7, §IV]. Indeed, the ring of invariants for plane cubics is 
generated by two polynomials of degrees 4 and 6, and the condition for a plane cubic to 
have a cusp is that both of these invariants vanish [17, §10.3]. This implies that, given seven 
general points in P 2 , there exist precisely 24 representations (4.4) of that configuration up to 
projective transformations. The net of cubics through these seven points defines a morphism 
from the corresponding del Pezzo surface of degree 2 onto P 2 . The branch locus is a quartic 
curve, and the 24 cuspidal cubics in that net correspond to the 24 flexes of that quartic. □ 

Remark 4.3. The coordinate system d and the connection to cuspidal cubics are due to 
Bramble [7]. Bramble also gives an explicit formula for the plane quartic in terms of W{E 7 )- 
invariant polynomials in di, . . . , d 7 [7, §VI] (there the coordinates are called t±, . . . , t 7 ). We 
first learned about the matrix D from [30, §6]. □ 
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Remark 4.4. Here is another way to see that the degree of 7 is 24. Form the 2 x 15-matrix 



The first row consists of our degree 7 polynomials in ci, . . . , C7, and the second row is the 
image under 7 of a random point C G P 6 . The corresponding fiber of 7 is defined by the 
ideal of 2 x 2-minors. One needs to verify that its scheme is reduced and consists of 24 
points in P 6 . In practice, this is done as follows. With probability one, we have R 7^ 0, 
and we can consider the affine scheme defined by r = 1, s oi = S 00 i/R, . . . ,t m = T m /i?. 
Using Macaulay 2 we verify that this is a 0-dimensional scheme of degree 168. This gives 
the number of solutions of (ci, . . . ,cj) G A 7 that belong to the fiber of this point. But if 
Aci, . . . , \cj were also a solution, we would need A 7 = 1 because all of the functions used to 
define 7 are of degree 7. Hence the degree of the fiber of 7 is 168/7 = 24, as desired. □ 

Remark 4.5. As indicated in the proof of Theorem 4.1, the map 7 is closely related to a 
map from the moduli space of 7 points in P 2 to the moduli space of plane quartics. Here, 
the reflection hyperplanes correspond to the 63 discriminantal conditions which give singular 
plane quartics. It follows that the non-hyperelliptic locus of the moduli space .4.3(2) is the 
image in the Gopel variety Q = .4.3(2) of the complement of the 63 reflection hyperplanes in 
P 6 . The fiber of size 24 corresponds to a choice of flex point on the plane quartic. 

In Section 6 we shall define an embedding of the Gopel variety Q into a certain 35- 
dimensional toric variety T . It lives in P 134 and H^(E7) acts on this embedding by permuting 
coordinates. That toric variety T and its polytope elucidate the combinatorial structure of 
the compactification Q = ^3 (2). The image of the 63 reflection hyperplanes give the 63 
boundary divisors in Q \ ^.3(2), and give 63 distinguished facets of the polytope of T. □ 

After completing the present work, we learned that our parametrization of the Gopel 
variety, as well as some of the material in Section 6, had already been found by Colombo, van 
Geemen and Looijenga in [12]. In particular, the indeterminacy locus of our map 7: P 6 ---> Q 
is determined in [12, Proposition 4.18]. See also [12, §1, Quartic curves] for an interpretation 
of this result in terms of the geometric invariant theory of plane quartics. We now state their 
result, and we strengthen it by including the answer to [12, Question 4.19]. This underlines 
the benefit of combining theoretical studies of moduli spaces with hands-on computations. 

Recall that a flat of a hyperplane arrangement is the intersection of some subset of the 
hyperplanes. The dimension of a flat will refer to the dimension of its projectivization. 

Theorem 4.6. (cf. [12]) The indeterminacy locus of 7 is the reduced union of 315 two- 
dimensional flats and 336 one- dimensional flats of the reflection arrangement of E 7 in P 6 . 
They correspond to the root subsystems of type D 4 and A 5 listed in row 9 and row 15 of 
Table 2 in Section 9. The ideal (r, is the intersection of these 651 linear ideals. 

Proof. The set-theoretic version of this theorem was proved in [12, Proposition 4.18]. We 
established the ideal-theoretic statement with a Macaulay 2 computation that is posted in 
our supplementary files. The following remark offers a more detailed explanation. □ 
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Remark 4.7. The indeterminacy locus of 7 is the zero set in P 6 of the ideal (r, s,,£.) = 
(r, sooi, . . . , s m , tooi, • • • , *m) hi the polynomial ring Q[c u c 2 , . . . , c 7 ] = Q[di, d 3 , . . . , d 7 ]. The- 
orem 4.6 states that this ideal is radical, and that it is the intersection of 651 = 315 + 336 
ideals which are generated by linear forms. One of these components is the height 4 prime 

{di—d 2 , dx—d 3 , di+d 4 +d 5 , d 2 +d 4 +d b , d 3 +d 4 +d 5 , di+d§+d 7 , d 2 +d e +d 7 , d 3 +d G +d 7 , , , 
d 2 -d 3 , d 2 +d 3 +d 4 +d 5 +d 6 +d 7 , di+d 3 +d 4 +d 5 +d 6 +d 7 , di+d 2 +d 4 +d rj +d 6 +d 7 ) . 

The 12 listed generators form a root subsystem of type D 4 in the root system E 7 . There are 
315 such subsystems in E 7 . By [39, §4], they are in bijection with the syzygetic triples of 
Steiner complexes, and also with the isotropic planes in (F 2 ) 6 . The latter will be discussed in 
Section 6, and a census of all root subsystems will be given in Table 2. Each root subsystem 
D 4 consists of 12 of the 63 reflection hyperplanes for E 7 , and these intersect in a P 2 inside 
P 6 = P(h). The rank 4 matroid given by (4.5) is the Reye configuration in [39, Figure 2]. 

The other 336 components of (r, s„t m ) correspond to an orbit of root subsystems of type 
A 5 . Each such prime ideal is linear of height 5, and it contains 15 of the 63 linear forms. 

To verify the claim, we intersect all 651 linear ideals in Macaulay 2. We note that a 
naive intersection may run for several hours and still not terminate, but a careful choice of 
ordering will finish in a matter of seconds. First, we use the c-coordinates rather than the 
d-coordinates. Each of our ideals contains at least one of the variables q. We separate them 
into seven groups of 93 based on which q they contain (there is some redundancy). Then 
we intersect each of these groups. Finally, we intersect the seven resulting ideals, and the 
result is (r, s,,t 9 ). We also used this procedure to get the ideals of the D4 flats and the A 5 
flats. We list the graded Betti tables for the indeterminacy locus of 7, the reduced union of 
the 315 flats of type D 4 , and the reduced union of the 336 flats of type A 5 , respectively: 



012345 012345 012345 

total: 1 15 84 168 126 28 total: 1 85 210 315 210 21 total: 1 36 315 595 420 105 

0: 1 0: 1 0: 1 

1: 

6: . 15 ... . 6: . 15 ... . 2: 

7: 7: 3: 

8: 8: . 70 210 . . . 4: . 21 ... 

9: . . 84 . . . 9: ... 315 210 . 5: 

10: ... 168 105 21 10: 21 6: . 15 315 595 420 105 

11: 

12: . . . . 21 . 
13: 7 



At this point, it is important to note that (r,s,,t,) has an alternative generating set, 
consisting of 135 products of linear forms defining reflection hyperplanes. We shall present 
them in Section 6. This fact makes it obvious that the variety consists of flats of E 7 . Colombo 
et al. used these 135 generators to prove the set-theoretic result in [12, Prop. 4.18]. □ 

5 Equations defining the Gopel variety 

In this section we first determine the equations and graded Betti numbers of the Gopel 
variety Q. Subsequently, we compute the prime ideal of the universal Coble quartic over Q. 
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This is the 12-dimensional subvariety in P 14 x P 7 whose fibers over Q are the quartics (2.12). 



Theorem 5.1. The six- dimensional Gopel variety Q has degree 175 in P 14 . The homogeneous 
coordinate ring of Q is Gorenstein, it has the Hilbert series 



1 + 8z + 36z 2 + 85z 3 + 36z 4 + 8z 5 + z e 



(5-1) 



and its defining prime ideal is minimally generated by 35 cubics and 35 quartics. The graded 
Betti table of this ideal in the polynomial ring Q[r, sqoi, • • • , ^m] in 15 variables equals 



8 



012 3 4 567 
1 70 609 1715 2350 1715 609 70 1 
1 

35 21 

35 588 1715 2350 1715 588 35 

. 21 35 



total 

1 
2 
3 
4 
5 
6 

Remark 5.2. Before giving the proof, let us point out some geometric consequences of this 
theorem. Since the degree of the numerator is less than the degree of the denominator in 
the Hilbert series, we see that the Hilbert function agrees with the Hilbert polynomial for 
all nonnegative inputs. One can calculate directly that the Hilbert polynomial of Q is 



35 fi 35 = 287 4 133 , 343 2 

-t 6 + — t 5 + 1 4 + 1 3 + 1 2 

144 48 72 
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-t+l. 



144 48 144 48 72 2 
We see from the self-duality of the Betti table that the canonical bundle of Q is 



UJg 



£xt S (Og,OM-^)) = Og(-l). 



Hence Q is a Fano variety. Since Q is arithmetically Cohen-Macaulay, and the Hilbert 
polynomial and the Hilbert function agree, this implies H*((?, Og(d)) = if either d > and 
i > or if d < and i < 6 [20, Cor. Al.15]. □ 

Proof of Theorem 5.1. Our point of departure for the derivation of Theorem 5.1 is [18, §IX.5, 
Proposition 8] and the subsequent corollary which gives a list of 63 cubics. We now review 
this construction. For each of the 63 non-zero half-periods e/2 G A/2A, we consider the 
linear involution on P 7 that is induced from the action on second order theta functions seen 
in (2.10). The fixed point set of this involution on P 7 is the union of two three-dimensional 
planes and H~, and each of these planes intersects the Coble quartic in a Kummer 
surface. The relations (2.14) on the coefficients of these two Kummer surfaces give the same 
cubic relation on r, Sijk,Ujk, and this cubic then lies in the ideal of the Gopel variety Q. 

For a concrete example consider e = (1,0,0). The corresponding involution (2.10) fixes 
the coordinates x oo ; ^ooi; ^oicb x ou an d it switches the sign on x\qq, Xi i, X\ W , xm. The 3- 
planes Hf are obtained by setting either of these two groups of four variables to zero. This 
specialization erases ten of the 15 summands in the formula for F T in (2.12). What remains 
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are the terms with coefficients r, s 001 , s 010 , Son^ioo- This represents a quartic Kummer sur- 
face in H+ resp. in H~ . The relation (2.14) holds for that Kummer surface. The argument 
in the previous paragraph shows that the following cubic lies in the ideal of Q: 

16r 3 + rt\ m + 4s 00 iSoioSoii - 4rs 2 )01 - 4rs 2 110 - 4rsg u . 

For a more complicated example consider the half-period e = r • (1, 0, 0)*. The three-plane 
H~ is defined by Xooo - ^100 = ^ooi — ^101 = ^oio - x iio — x on ~ x i\i — 0, and the three-plane 
is defined by xooo + ^ioo = ^ooi + x wi = ^oio + £110 = £011 + = 0. The cubic relation 
on the Kummer surface obtained by restricting F to H~ or equals 

16(s ii + 2r) 3 + (son + 2r)(2i uo + 2t 10 i + 2t 00 i + 2t io) 2 
+4(2s 00 i + 2s io + tioo)(2sioo + *on + 2s in )(2s 101 + t ni + 2s no ) 
-4(s ii + 2r)(2s 00 i + 2s io + ^ioo) 2 - 4(s n + 2r)(2s 100 + t 011 + 2s m ) 2 

-4(s ii + 2r)(2s 101 + t in + 2s no ) 2 . 

In this manner we obtain 63 cubic equations for the Gopel variety Q, as stated in [18, 
Corollary on p. 186]. However, only 35 of these are linearly independent. A vector space 
basis consists of those cubics that come from the 35 nonzero even theta characteristics. 

In addition to the 35 cubics constructed from Kummer surfaces as above, the ideal of Q 
has 35 minimal generators of degree 4. Here is an example of such a quartic generator: 

48r 2 sioiSno — 12sq 11 s 10 iSho — 12s 2 00 sioiSno — 4s 2 )01 sioo ,s m — 4sq 10 SiooSih + 8s 2 )11 Sioo s iii+ 
4sio siii-16rsooiSioiSm + 8s ioSoiiSioiSiii + 8siooSioi s iii-l 6rs oioSiioSiii+8s oiSoiiSiioSiii 
+8sioos no sm + SsooiSoio^in— 16rs ns ni — 4sioiSnos ni — 4siooSin— sioo-smtooi — s ioo s m^oio 
— 2s oi s oio^on + 4rson^on + s ioi s iio^oii + ^looSiii^oii+^ioo^iii^ioo+^ioo^iii^ioi + s ioo s iii^iio 

— Sioo^on^ioi^no + s iio^oio^ioi^in + sioi^ooAio^m + 2sooi s oio^in — 4rs nt ni — 2s 10 iSnotin 

The explanation for the derivation of such quartics will come in the proof of Theorem 6.1. 
From the 70 minimal generators, a Grobner basis for the ideal of Q can be computed in 
Macaulay 2, and from this one finds the degree 175 and the Hilbert series in Theorem 5.1. 

The group Sp 6 (F 2 ) has exactly two distinct irreducible representations of dimension 35, 
and these two representations are given respectively by the cubic generators and the quartic 
generators of our ideal Q. In particular, we can obtain all 70 generators by lifting the two 6x6- 
matrices /x and v in Section 4 to Q[r, S001, • • • , £111] and applying these to one representative 
cubic and one representative quartic, for instance those displayed above. In Section 6 we shall 
present an alternative model for these two Sp 6 (F 2 )-submodules Q[r, s oi, • • • , tin], namely in 
terms of binomials in a polynomial ring in 135 unknowns. 

It remains to be proved that the 35 cubics and the 35 quartics actually generate the prime 
ideal of Q. This is what we shall explain next. The 70 generators, as well as explicit matrix 
representations for /j, and v on the 15 Coble coefficients, are available in our supplementary 
materials, and the reader can use these to verify all of our calculations in Macaulay 2. 

The main idea is to show that the scheme defined by this ideal is generically reduced 
and satisfies Serre's criterion (Si). These conditions ensure that our ideal is radical by [19, 
Exercise 11.10]. By the earlier results of Coble [11, §49] and Dolgachev-Ortland [18, §IX.7, 
Theorem 5], we already know that the cubics alone cut out Q set-theoretically, and hence 
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the variety of our ideal is irreducible of codimension eight. So, to show 'generically reduced' 
only requires showing that at some point on Q, the Jacobian matrix has rank 8. This is 
easily done by substituting randomly chosen points from P 6 into the map 7. 

To prove Serre's condition (Si), we show that our ideal defines an arithmetically Cohen- 
Macaulay scheme. This is done by adding seven random linear forms to the ideal generated by 
the cubics and quartics. The resulting ideal defines an artinian quotient of Q[r, s oi, • • • , ^m], 
and we verified that its Hilbert series is the numerator seen in (5.1). The computation shows 
that the seven linear forms are a regular sequence modulo our ideal, and hence we get that 
the scheme defined by the 35 cubics and 35 quartics is arithmetically Cohen-Macaulay. 
We remark that the regular sequence computation is intensive and took approximately 15 
minutes to perform. An alternative method is to work over the finite field Z/101. The Hilbert 
series remains unchanged, but the regular sequence computation takes only a few seconds to 
perform. This implies the result via standard flatness and semicontinuity arguments. 

At this point we have shown that the ideal generated by our cubics and quartics is prime. 
To infer the Gorenstein property, we use Stanley's result [49, Theorem 4.4], which states that 
a Cohen-Macaulay domain whose Hilbert series has a palindromic numerator is Gorenstein. 

We finally derive the Betti table asserted in Theorem 5.1. We have already shown that the 
Gopel variety Q is arithmetically Gorenstein, which implies that the Betti table is symmetric 
[19, Corollary 21.16]. Multiplying the numerator of the Hilbert series (5.1) by (1 — z) 8 gives 
the graded Euler characteristic (or K-polynomial) of the minimal free resolution of Q. By 
the Cohen-Macaulay property, the Castelnuovo-Mumford regularity of Q is the degree of 
the numerator of the Hilbert series, which is 6. Therefore, we have the following partial 
information about the Betti numbers: 
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: 1 


70 
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70 
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e 


b 








35 
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21 
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where the unknown entries satisfy the linear equations 

a- c = 588,b-d + f = -1715, -2e + g = 2350. 

Now, there are 21 linear syzygies on the 35 cubic generators, and we compute the corre- 
sponding 35 x 21-matrix whose entries are linear forms in Q[r, sooi, • • • , £m]- Evaluating this 
matrix at a random point in Q 15 shows that it has full rank 21 over the rational function 
field Q(r, sooi, • • • , im)- Hence its columns are linearly independent. This implies c = 0, and 
consequently f = e = b = 0. The above linear equations then imply a = 588, d = 1715, and 
g = 2350, so the Betti table is as claimed. This concludes our proof of Theorem 5.1. □ 

Now we study the universal Coble quartic C in Q x P 7 . This twelve- dimensional variety 
is defined as the closure in P 14 x P 7 of the set of pairs ((r : s oi : • ■ ■ : tm), (xooo : • • ■ : £111)) 
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such that (r : s oi : ■ • • : t in ) is a point in the non-hyperelliptic locus of .4.3(2) C Q, and the 
pair is a point on the hypersurface of bidegree (1, 4) that is given by the polynomial (2.12). 

Corollary 5.3. The bi- homogeneous prime ideal of the universal Coble quartic in P 14 x P 7 is 
minimally generated by 71 polynomials, namely the 35 + 35 equations of bidegrees (3, 0) and 
(4,0) for the Gopel variety as in Theorem 5.1, along with the bidegree (1,4) equation (2.12). 

Proof. The equation (2.12) cuts out a codimension one subscheme C of Q x P 7 . Since Q is 
arithmetically Cohen-Macaulay in its embedding in P , this implies that (2.12) is a non- 
zerodivisor on Q x P 7 , and that the scheme C is arithmetically Cohen-Macaulay as well. To 
show that our ideal is prime, it suffices to show that it is irreducible and generically reduced 
(as explained in the proof of Theorem 5.1). This will imply C = C. 

Consider the projection map tt: C — > Q. The fibers of it have constant dimension 6 since 
the Coble coefficients cannot simultaneously vanish on any point of Q. Let IA C .4.3(2) denote 
the non-hyperelliptic locus. Over the subvariety U C Q, the fibers of 7r are Coble quartics, 
and hence irreducible. Since C is arithmetically Cohen-Macaulay, all of its irreducible 
components have the same dimension, namely 12. The closure of n^iU) gives an irreducible 
component. The complement of tt~ 1 (U) in C is the preimage of QXU. But dim7r -1 (C7\£V) < 
11, so it cannot be an irreducible component. Hence the variety C is irreducible. 

Now consider the Jacobian matrix of the 71 equations defining C. Since the 70 equations 
for Q involve no x coordinate, this matrix can be put in block triangular form. Choose any 
point ((r : s oi : • • ■ : *m), (^000 : • ■ ■ : £111)), where (r : s oi : ■ • • : t in ) G A 3 (2) is a non- 
singular point of Q, and (xooo : • ■ • : £111) is a non-singular point on its Coble quartic. Then 
this is a non-singular point of C . Hence the scheme C defined by the 71 proposed equations 
is also generically reduced. We conclude that they generate the prime ideal of C = C . □ 

6 A toric variety for seven points in P 2 

In this section we present an alternative embedding of the Gopel variety, with beautiful com- 
binatorics. We learned this from Dolgachev and Ortland [18, §IX.7]. Here, the Gopel variety 
sits in the high-dimensional projective space P 134 whose coordinates are the 135 Gopel func- 
tions. The corresponding ideal is generated by linear trinomials, cubic binomials and quartic 
binomials. We will first describe the combinatorial structure of the Gopel functions. A future 
project will be devoted to the analogous toric varieties for other Macdonald representations. 

Consider the six-dimensional vector space (F 2 ) 6 over the two-element field F 2 . We fix the 
following non-degenerate symplectic form on this 64-element vector space: 

(x,y) = x ± y 4 + x 2 y 5 + x 3 y 6 + x^y x + x 5 y 2 + x 6 y 3 . (6.1) 

A linear subspace V C (F 2 ) 6 is isotropic if (x, y) — for all x, y e V. Clearly, all subspaces of 
dimension < 1 are isotropic, and there are no isotropic subspaces of dimension > 4. There are 
precisely 315 isotropic subspaces of dimension two, and 135 isotropic subspaces of dimension 
three. The former are called isotropic planes, and the latter are called Lagrangians. There 
are 63 non-zero vectors in (F 2 ) 6 , and each Lagrangian contains seven of these. Each isotropic 
plane contains three of these, and is contained in precisely three Lagrangians. 
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Table 1: Cayley's bijection between lines in (F2) 6 and positive roots of E7. 

The root system E 7 consists of 126 vectors in a seven-dimensional inner product space. 
If we take this space to be the hyperplane in M 8 given by coordinate sum zero, then the 
63 = (2) + (3) positive roots are — ej for 1 < i < j < 8 and |(e 8 + ^ igo .ej — X^cr e j) 
where a runs over all three-element subsets of {1, 2, ... , 7}. We label the 63 = (2) + ( 3 ) 
positive roots of E 7 by the pairs in {1,2,..., 8} and the triplets in {1,2,..., 7}. 

We copied Table 1 from Cayley's paper [8]. It fixes a particular bijection between the 
positive roots of E7 and the vectors in (F 2 ) 6 \{0}. Its rows represent the first three co- 
ordinates and its columns represent the last three coordinates of a vector in (F2) 6 . For 
instance, the triple 247 corresponds to the vector (0, 0, 1, 1, 1, 0), and the pair 24 corresponds 
to (1,1,1,1,0,0). This bijection has the property in [18, Lemma IX. 8]: two positive roots 
in E7 are perpendicular in M 8 if and only if the corresponding vectors x,y G (F 2 ) 6 satisfy 
(x, y) = in F 2 . Cayley had constructed his table in the 1870s to have precisely this property. 

The orthogonality preserving bijection between (F 2 ) 6 \{0} and the positive roots of E 7 
shows that the Weyl group of E 7 (modulo its two-element center) is isomorphic to the group 
Sp 6 (F 2 ) of 6 x 6-invertible matrices over F 2 that preserve the symplectic form (6.1). The 
order of that symplectic group is easily seen (cf. [11, §11.23, (7)]) to be 

(2 6 - 1) • 2 5 • (2 4 - 1) • 2 3 • (2 2 - 1) • 2 1 = 36-8! = 1451520. (6.2) 

With the relabeling given by Cayley's table, the 135 Lagrangians fall into two classes with 
respect to permutations of the set {1, 2, 3, 4, 5, 6, 7}. First, there are 30 Fano configurations 

{124, 235, 346, 457, 561, 672, 713}. (6.3) 

We denote this configuration by /i 2 34567- Second, there are 105 Pascal configurations like 

{ 12, 34, 56, 78, 127, 347, 567}. (6.4) 

We denote this configuration by c/7123456 • If the configurations are changed by a permutation 
of {1,2,3,4,5,6,7} then the indices of the label /. or g, are permuted accordingly. Hence 
the labeling is not unique. For example, the Fano configuration /i 2 34567 in (6.3) can also 
be labeled / 2 34567i, and the Pascal configuration #7^3456 in (6.4) can also be labeled #7345612- 
Some permutations give rise to a sign change in /. or g,. 
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The 315 isotropic planes in (F 2 ) 6 fall into three classes. Each class has cardinality 105. 
Using the labeling given by Cayley's Table 1, representatives of these three classes are 

{12, 123, 38} and {123, 145, 167} and {12, 34, 567}. (6.5) 

Consider the six-dimensional variety of unlabeled configurations of seven points in the 
projective plane P 2 . There are natural correspondences, described in [18, §IX], that take 
such configurations to Cayley octads, hence to plane quartic curves, and hence to abelian 
threefolds. Here is how to make this completely explicit. We associate with each Lagrangian 
in (F2) 6 a homogeneous symmetric polynomial in the brackets [ijk], which are the 3x3- 
minors of the 3 x 7-matrix of homogeneous coordinates on (P 2 ) 7 . These bracket polynomials 
are called Gopel functions in [11, 18]. Each of the 30 Fano configurations translates verbatim 
into a bracket monomial. For instance, the Fano configuration (6.3) translates into 

/i234567 = [124] [235] [346] [457] [561] [672] [713] . (6.6) 

The 105 Pascal configurations represent six points lying on a conic, e.g. let Q7 denote the 
quartic bracket polynomial that vanishes when the points 1, 2, 3, 4, 5, 6 lie on a conic: 

Q 7 = [134] [156] [235] [246] - [135] [146] [234] [256]. (6.7) 

Then the Pascal configuration (6.4) translates into the bracket binomial 

37123456 = [127][347][567]-Q 7 . (6.8) 

We have thus defined a coordinate system on P 134 , consisting of 30 coordinates /. and 105 
coordinates <?., and we have defined a rational map 

( P 2)7 pl34 ( 6 _ g) 

whose coordinates are multi-homogeneous bracket polynomials of degree (3,3,3,3,3,3,3). 
This map factors through the six-dimensional space of configurations of seven points in P 2 . 

One desirable property of this coordinate system is that it fits well with E 7 . Consider the 
action of the symplectic group Sp 6 (F 2 ) on the 135 Lagrangians in (F 2 ) 6 . Then the induced 
action on our 135 Gopel coordinates /, and g, is realized by signed permutations. 

Dolgachev and Ortland [18, §IX.7, Proposition 9] describe the following 315 linear re- 
lations among the Gopel functions. For each of the 315 isotropic planes, there is a linear 
relation among the three Gopel coordinates indexed by the three Lagrangians containing 
that plane. For instance, the three isotropic planes in (6.5) determine the linear relations 

#3124567 — #3124657 + #3124756 #1234567 — /l243765 + /l243675 #5123467 _ #6123457 + #7123456 



{12, 38, 123, 45, 67, 345, 367} {123, 145, 167, 18, 23, 45, 67} {12, 34, 567, 58, 67, 125, 345} 
{12, 38, 123, 46, 57, 346, 357} {123,145,167,247,256,346,357} {12, 34, 567, 68, 57, 126, 346} 
{12, 38, 123, 47, 56, 347, 356} {123,145,167,246,257,347,356} {12, 34, 567, 78, 56, 127, 347} 

In each column we list the corresponding triple of Lagrangians. These linear forms arise from 
the three-term quadratic Pliicker relations among the brackets [ijk]. These linear forms span 
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a 120- dimensional space, so they cut out a 14- dimensional linear subvariety L C P 134 . Thus, 
the rational map (6.9) factors through this linear space: 

(P 2 ) 7 L C P 134 . (6.10) 

Now, returning to the Coble quartic F T in (2.12), its 15 coefficients span the same irre- 
ducible Sp 6 (F 2 )-module as the 135 Gopel functions. Coble states in [11, §IV.49,(16)] that 

r ; s 00l5 s 01Cb s 01l5 s 10Cb Sioi , SjlO) s lll> ^001) ^010) ^011)^100) ^101: ^110: ^111 

can be expressed as integer linear combinations of the Gopel functions /. and g,. 

Constructing such linear combinations turned out to be a non-trivial undertaking, but 
we succeeded and found an explicit list of 135 transformation formulas such as 

/i234657 = — 4r — 2sooi — 2s io — 2son — 2sioo ~~ 2sioi — 2sno — 2s m 

— ^001 — ^010 — ^011 ~~ ^100 ~~ ^101 — ^110 — ^111 /g -q\ 

/l237654 = — 4sioi — 8r — 4s oi ~~ 4Sioo ~~ 2t io 

Of course, these 135 linear forms in the 15 coefficients of F satisfy the 120 linear trinomials 
above. We also have explicit formulas for writing r, s, t in the Gopel coordinates g m . The 
complete list of transformation formulas is posted on our supplementary materials website. 

It is important to note that the change of basis in (6.11) is unique up to scalar multi- 
ple, since we require it to be Sp 6 (F 2 )-equivariant, and the 15- dimensional Sp 6 (F 2 )-module 
appears with multiplicity 1 in the permutation representation given by the Gopel functions. 
Consequently, the expressions in (6.11) are not just some random transition maps. In finding 
them explicitly, we resolved an issue that was left open by Coble in [11, §IV.49, (16)]. 

Now we are ready to connect the map (6.9) with the map (4.1). 

Theorem 6.1. The closure of the image of the map (P 2 ) 7 P 134 described prior to (6.9) 
lies in a 14- dimensional linear subspace of P 134 and it is linearly isomorphic to the Gopel 
variety Q in P 14 . Its prime ideal in the polynomial ring in 135 unknowns f, and g, is 
minimally generated by 120 linear trinomials, 35 cubic binomials, and 35 quartic binomials. 

Proof. The transformation (6.11) defines an isomorphism L ~ P 14 that is compatible with 
the rational map (P 2 ) 7 —• » L in (6.10) and the rational map (P 2 ) 7 — » Q C P 14 in (4.1). 
Therefore, the image of our map (P 2 ) 7 — -> P 134 in (6.9) is linearly isomorphic to the Gopel 
variety Q. Moreover, under this projective transformation, the ideal of the image of the 
rational map (P 2 ) 7 — -» P 134 in L is mapped to the ideal of Q in P 14 , which is the Gorenstein 
prime ideal generated by the 35 cubics and 35 quartics in Theorem 5.1. 

It remains to be shown that the cubic and quartic generators can be represented by 
binomials modulo the 120 trinomials. The existence of binomial relations comes from the 
fact that Q sits in the image of a monomial map, to be defined below. Using that, we will 
now give binomial relations, whose existence is explained in Theorem 6.2. There are 630 
cubic binomials of the form 

5 , 51324675 l 6123457/l235746 — /l236745fl , 5123467fi , 6132457- (6.12) 
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They form a single orbit under the Sp 6 (F 2 )-action by permuting Gopel coordinates. We can 
select a suitable subset of 35 of them to serve as generators for the ideal of Q C P 134 . 

Further, there are precisely 12285 quartic binomials that are not in the ideal generated 
by the 630 cubics. These form two Sp 6 (F2)-orbits of sizes 945 and 11340. Representatives 
for these two orbits are, respectively, 



/l234567,/l234675#1243657#2143756 — /l234576,/l234765#1243756#2143657) 
/l234576 /l234765 /l235647 /l236457 — /l234675 /l234756 /l235467 /l236547 • 



(6.13) 



We select a suitable subset of 35 of these quartics as well. So, now we have 70 binomials 
that vanish on the image of (6.10). A computation verifies that the image under (6.11) of 
the ideal they generate coincides with the ideal of Q C IP 14 in Theorem 5.1. □ 

We now present a completely explicit formula for the map P 6 Q C P 134 . It factors 
through the configuration space of seven points in P 2 as follows. A point {d\ : d 2 '■ ■ ■ ■ '■ d 7 ) in 
P 6 is sent to the configuration given by the columns of the matrix D in (4.4). We specialize 
the 35 Plucker coordinates [ijk] to the 3 x 3-minors of D. The result is the product 

[ijk] = (di - dj)(di - d k )(dj - d k )(di + dj + d k ). 

Under this specialization, the quartic bracket polynomial Qi for the condition (6.7) that six 
points lie on a conic maps to a product of 16 linear forms. Substituting these expressions 
into the Gopel functions (6.6) and (6.8), we obtain a list of 135 polynomials in di, d 2 , • • • , d 7 , 
each of which is a product of 28 linear forms. All 135 such products share a common factor 
of degree 21, namely, n^-C^ ~ dj). Removing that common factor, we recover the formulas 
for the Gopel coordinates as products of seven linear forms. For instance, we find 

/i234567 = (di+d 2 +d 4 ) (d 2 +d 3 +d 5 )(d 3 +d 4 +d 6 ) (d 4 +d 5 +d 7 ) (d 5 +d 6 +di)(d 6 +d 7 +d 2 )(d 7 +d 1 +d 3 ) , 
#7123456 = (di+d 2 +d 7 ) (d 3 +d 4 +d 7 ) (d 5 +d 6 +d 7 ) (d 1 -d 2 )(d 3 -d 4 ) (d 5 -d 6 )(-di-d 2 -d 3 -d4-d 5 -de) 

If we now substitute these products of linear forms into the linear map that inverts (6.11) 
then we recover the map P 6 P 14 from (4.1). These formulas also appear in [12, §3]. 
We now introduce variables for the 63 reflection hyperplanes of E 7 . Denote 



J. ; 



di - (di + ■ ■ 
di + dj + d k . 



+ d 7 ), 



(6.14) 



Then, we obtain 



/l234567 = ^124^137^156^235^267^346^457) /l234576 
/l234657 — ^124^137^156^236^257^345^467) /l234675 

fl 243765 

#1234567 = x I x 2 / 3X^Xq 7 Xi 23 X\^Xiq 7 , #1234657 
#1234756 = ^1^23^47^56^123^147^156) #1243567 

#7162534 



^124^136^157^235^267^347^456) 
X iuX 135^167X236^257^347^456 j 

x l x 2 l i x AQ x b7 x l2Z x lA& x lb7-i 
2:1X24X35X67X124X135X167) 
X7Xi6X25X34Xi67X257X347- 



(6.15) 
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Each of these x-monomials is squarefree of degree 7, and represents one of the 30 Fano 
configurations like (6.3), or one of the 105 Pascal configurations like (6.4). 

This is the moment when we come to the object promised in the title of this section. If 
we regard Xi, Xy, Xijk as formal variables, the formulas in (6.15) define a monomial map 



The closure of the image of m is a toric variety T in P 134 . We call T the Gopel toric variety. 

Let A denote the 63 x 135-matrix with entries in {0, 1} representing the monomial map 
m. The rows of A are labeled by the 63 parameters Xi,Xij,Xijk, or, using Cayley's bijection, 
by the 63 vectors in (F 2 ) 6 \{0}. Here we erase all occurrences of the index "8" in Table 1 
so as to identify (F2) 6 \{0} with the subsets of {1, 2, 3, 4, 5, 6, 7} that have cardinality 1, 2 or 
3. The columns of the matrix A are labeled by the 135 Gopel coordinates f,,g,, and thus 
by the 135 Lagrangians in (F2) 6 . Each row of A has 15 entries 1; the other 120 entries are 
0. Each column of A has 7 entries 1; the other 56 entries are 0. By computing the Smith 
normal form, one checks that our 63 x 135-matrix A has rank 36 over any field. 

The convex hull of the columns of A is a convex polytope of dimension 35 in IR 63 . We 
denote this polytope also by A and we call it the Gopel polytope. The Gopel polytope A has 
135 vertices and 63 distinguished facets, each given by the vanishing of one of the coordinates 
in the ambient IR 63 . However, these are not all of the facets of A. In fact, it is a difficult 
computational problem, which we could not solve, to determine the number of facets of 
A. Each distinguished facet contains 135 — 15 = 120 of the vertices, and is indexed by a 
vector v € (F2) 6 \{0}. Thus each distinguished facet of A corresponds combinatorially to the 
complement of a Lagrangian in (F2) 6 \{0}. These facets are indexed by the 63 x- variables. 

Theorem 6.2. The Gopel toric variety T has dimension 35 in P 134 . It has 63 distinguished 
boundary divisors, each given by the vanishing of one parameter The Gopel 

variety Q is the (ideal-theoretic) intersection of T with the linear space L ~ P 14 in (6.10). 

Proof. By construction, the Gopel polytope A is the polytope associated with the projective 
toric variety T, and so the first assertion follows from the above discussion. The ideal for 
the embedding T C P 134 is the toric ideal associated with the matrix A. We computed all 
minimal generators of this ideal of degree at most 4. There are no linear forms or quadrics 
in the toric ideal of A, but there are 630 cubic binomials (6.12) and 12285 quartic binomials 
(6.13) that are not in the ideal generated by the 630 cubics, which are straightforward to find. 
We have already shown that the cubics and quartics generate the ideal of Q inside L. □ 

It would be desirable to find a Markov basis of T, i.e. a set of minimal generators of the 
toric ideal 1^ of the matrix A. This turns out to be nontrivial. In fact, we need binomials 
of up to degree at least 6 to generate the toric ideal. We constructed examples of binomials 
of degree 5 and 6 that lie in but are not generated by elements of lower degrees: 

/l2345765 , 21637455 , 41326575 , 61425375 , 7132456 /l2365745 , 21437565 , 41625375 , 61324575 l 7132645 
/l234576/l2436575 , 21534675 , 31527465 , 41625375 l 6123547 — Zl236475,/l2435765 , 21537465 , 31625475 , 4123567fi , 6152734 

Both of these binomials are indispensable: they represent two-element fibers of the semigroup 
map A: N 135 —> N 63 . This implies that they must appear in every Markov basis of I_a. 
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We also do not know whether T is projectively normal, or arithmetically Cohen-Macaulay. 
If this holds then T is the toric variety of the polytope A in the strict sense of [13, §2.3]. 
Also we do not know the Hilbert polynomial (or Ehrhart polynomial). At the present time, 
our 63 x 135-matrix A seems too big for the packages 4ti2 (for computing Markov bases), 
polymake (for computing facets), and normaliz (which is discussed in [13, Appendix B.3]). 

7 The universal Coble quartic in P 7 x P 7 

Recall that for each r G fiz corresponding to a smooth non-hyperelliptic curve of genus three, 
the Coble quartic C T is the unique hypersurface of degree 4 in P 7 whose singular locus is the 
Kummer threefold corresponding to that curve. The universal Coble quartic C(2,4) is the 
Zariski closure in P 7 x P 7 of the set of pairs (u, x) where u 6 $(4. 3 (2,4)) corresponds to a 
non-hyperelliptic curve and x lies in the Coble quartic hypersurface corresponding to u. In 
this section we derive generators for the ideal of C(2, 4) in the theta coordinates u and x. 

We shall derive formulas for the coefficients r, s 9 , t, of the quartic polynomial F T in (2.12) 
in terms of the second order theta constants u. Our formulas are fairly big, and they define 
an explicit rational map of degree 64 from the Satake hypersurface to the Gopel variety: 

s^lg. (7.i) 

This is realized concretely as a map P 7 P 14 by listing 15 polynomials r(u), s,(u), t.(u) 
in the variables u = (u 000 , • • • ? M m)- ft agrees with the abstract map (3.2) over the non- 
hyperelliptic locus of A4 3 (2,4) and collapses each of the 36 components of the hyperelliptic 
locus to a single point (see, for example, (8.9)). We also determine the zero set of these 
polynomials, which is the indeterminacy locus of (7.1) in the hypersurface S = $(.4.3(2,4)). 
The following theorem summarizes the results to be proved in this section: 

Theorem 7.1. (a) The coefficients r, of the Coble quartic can be expressed as polyno- 
mials of degree 28 in the eight theta constants u. The resulting polynomial F from (2.12) 
is the sum of 372060 monomials of bidegree (28,4) in (u, x). 

(b) The locus in the Satake hypersurface S that is cut out by the coefficients r, s,,t. equals 
the Torelh boundary Sing(5) = tf(.4 3 (2, 4)) \ $(A1 3 (2, 4)) (cf. Proposition 3.2). 

(c) The prime ideal of the universal Coble quartic equals (F,S), where S is the Satake 
polynomial of bidegree (16, 0). Hence C(2, 4) is a complete intersection of codimension 2 
in P 7 x P 7 , and its bidegree equals 16U(28U + 4X). Here we write the cohomology ring 
0/P 7 x P 7 as 1i[U, X]/(U 7 , X 7 ), with U and X representing the two hyperplane classes. 

Proof, (a) We will construct 15 polynomials of degree 28 in the eight unknowns uqoo, • • • , Mm- 
The polynomial r(u) has 5360 terms, each polynomial Sjjfc(u) has 7564 terms, and each 
polynomial iijk(u) has 8114. They are available on our supplementary materials website. 

We begin by outlining the derivation of these polynomials. Consider the Coble quartic 
C T corresponding to some non-hyperelliptic r. Let e be a non-zero length 3 binary vector, 
so e/2 G A T [2] represents a 2-torsion point on the abelian threefold A T . As in the proof of 
Theorem 5.1, we write Hf for the two three-dimensional fixed spaces of the involution (2.10) 
on P 7 defined by 02[cx](t; z) \-> O2 [ct](t; z + e/2). The fixed space H+ is given by 

x s = 0, S $ e 1 , (7.2) 
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where the orthogonal complements are taken with respect to the usual Euclidean form on 
(F 2 ) 3 . On the subspace ~ P 3 we use the homogeneous coordinates 



(^ooo : x 5l : x Ss : x Ss ), 



Si e 



Substituting (7.2) into F T , we get the equation of a Kummer surface 

s a ■ m a + t e ] [ x a = where m a = - ^ 



cnH: 



2 2 
X u X u+a- 



(7.3) 



<r6e J 



Its singular locus is H+ n Sing(C r ). We claim that one of the 16 singular points of is 



Pc 



4' 



(e 2 [000](r;-):e 2 [001](r; 



© 2 [lH](r; i )). (7.4) 



Since this point lies in the image of the Kummer map k t , it is by definition in the Kummer 
threefold K 3 = k t (A t ) = Sing(CV). That p e is in the fixed space can be seen from the 
transformation properties of the second order theta functions under A T [2]. Indeed, 



e 2 [5](r; 



9(2 



'5) ■ exp 



TT'l 
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Furthermore, 



9(2t;- + t5) 



for any 6 G (F 2 

r5)-e- (2t)S) 
- t5) - {l/2)8 t {2r)5) 



(7.5) 



9(2t;-C- + t5)) = 9(2r;C- 

= 6(2t;^ + t5) ■ exp 2vrz(5'(^ 

= 9(2t; I + t8) ■ exp [iritfe] . 

If 6 $ e\ then exp [7ri<5*e] = -I. So 9(2r; f +rS) = 0. Thus, p e is in H+ n K 3 . 

Now recall (from Example 1.1) that the equation of a Kummer surface can be expressed 
in terms of the coordinates of any of its 16 singular points. We apply this to the Kummer 
surface ■ Namely, let L\ denote the signed 4 x 4- minor of the corresponding matrix (1.1) 
which is obtained by deleting row 1 and column % + 1. Then the coefficients of (7.3) are 
proportional to Lq, L\, L\, L\, L\. Explicitly, we have the relations 
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Combining these formulas, we see that with respect to the normalization r = 1, the coeffi- 
cients of the Coble quartic (2.12) are given by 
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L 001 L ou L a 

s no = yfoT' s m = 7§iT ' and l " = for a G Ff\{0}. 

The are polynomials of degree twelve in the eight coordinates of the point p e in (7.4). 
The next step is to relate these coordinates G2[cx](r; |) to the variables u on our P 7 . 
We introduce the following notation for the first order theta constants: 

Tl, = 0[e|e'](r;O). (7.6) 

By the inverse addition formula (2.8), we have 

8& 2 [a](T; e -)e 2 [a + 5}(r; e -) = £ (-1)^7^, (7.7) 

e'e(Z/2)» 

In the special case 5 = 0, the left hand side becomes a square: 

8© 2 2 M(r;^) = £ (-1)^'^^, (7.8) 

e'G(Z/2)9 

Now, for < i < 3, the quantity L\ is actually a polynomial of degree 6 in 0|[*](t; f ). Hence 
formula (7.8) allows us to write these L\ as polynomials in the variables T°. For example, 

r 100 _ 1 /T^O T^O rpQ rpQ _ rpQ rpQ rpQ rpQ \ ,rpQ rpQ rpO rpO _ rpO rpO rpQ rpQ \ 

~~ Q4^ 000 010 100-* 110 1 001 1 011 1 101 1 111 A 1 000 J 001 -* 100-* 101 -*oio-*on-* 110 -* 111/ 
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rlOO _ 1 (rpQ rpO rpQ rpQ _ rpQ rpQ rpQ rpQ \frp0 rpQ rpQ rpQ , rpQ rpO rpQ rpQ \ 

~~ 32 00 10 00 001 011 101 111/ V-* 000 001 100 101 ' -*010 011 110 111/ 

,0 N 

110J' 

rlOO _ (rpO rpO rpO rpO _ rpO rpO rpO rpO \ ( rpO rpO rpO rpO _ rpO rpO rpO rpO \ 
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The quantity L\ can be written in the form 

L\ = HQ 2 [5](r; e -)-M^ (7.9) 

where M e is a polynomial of degree 4 in 0|[*] (r; |). Therefore, we can apply (7.8) to write M e 
as a polynomial of degree 8 in the T e e ,. To deal with expressions of the form n^e^ @2[^](' r ; f), 
we group them into two products of pairs of theta functions and apply (7.7). Note that some 
of the Tl, on the right hand side of (7.7), namely the ones such that e' G e x , correspond to 
odd characteristics and thus vanish. Hence we can write, for example, 



rlOO _ 1 ( TT rp0\ ( frp001\2(rp001\2 ( rpOOl \ 2 ^001 \ 2 \ 

4 _ o ' 1 11 J e'J ' ( l J 000 J l J 100 J _ l J 010 ) \ 1 no ) )■ 
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In this manner, we express all of the coefficients of the Coble quartic in terms of the T e e ,. 
Clearing the denominators, we get the following expressions of degree 28 for r, s a ,t a : 

r = n (n^-n^V ( 7 - iq ) 
*r=(-2)-(n i ?+ U T e)- n fn^-n^). ( 7 - n ) 

where W CT is a polynomial of degree 4 in the T e e , of the form {T*f{T*f - {T*f(T*f. 

To write these coefficients as polynomials in the u a , we'd like to use the addition formula 

(t:,) 2 = ^(-ir e v«. +e . (7.i3) 

In order to do this, we must write the coefficients of the Coble quartic as polynomials in 
the squares of the T e e ; . The expressions for r, s a ,t a above are not yet in this form. However, 
by expanding the products we observe that these expressions can be written as polynomials 
in the squares of the T e e , together with rLeF^ ^e°- Thus, it remains to express the quantity 
fLeF 3 as a polynomial in the squares of the T e e ,. To do this, we make use of the following 
special case of Riemann's theta relations [6, Exercise 7.9]: 

n t ' - n t ° = n ^ ^ 

Squaring this formula, we get 

n ra 2 + n ( r °) 2 - 2 n t ° = n ro 2 - 

e'6e ± e'^e ± e'€F| e'Ge ± 

Thus, we get a formula for n^eF 3 as a polynomial in the squares of the T variables. Then, 
we apply the formula (7.13) to the expressions (7.10), (7.11), (7.12) for r,s a ,t a . We obtain 
the 15 polynomials r(u), s CT (u), t CT (u) of degree 28 in the variables u that had been promised. 

(b) We next argue that, set-theoretically, the variety of the ideal (r,s a ,t a ) equals the 
Torelli boundary in the Satake hypersurface S. First, suppose that p £ S lies in the zero 
locus of (r,s a ,t a ). By (7.10) and Riemann's theta relations (7.14), the quantity 

e'Ge- 1 - e'ee x e'^e 1 - 

vanishes at p for at least one e £ \ {0}. By the proof of Proposition 3.2 or [26, Theorem 
3.1], in order to show that p lies in the Torelli boundary, it suffices to show that at least 2 
first order even theta constants T e e , vanish at p. We distinguish the following two cases: 
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Case 1. Suppose rie'ee- 1 - ^t 1 vanishes for two distinct e = ei,e2 G F| \ {0}. Since, these 
monomials have disjoint sets of variables, at least two first order even theta constants vanish. 

Case 2. Suppose that (Ile'ee- 1 - ^ — FL'^ vanishes for exactly one e 6 F| \ {0}. 
Then by (7.11), s e is a product of (H'e^ T ° + Ue'^ T °) with 

some nonzero factors. Thus, 
(Ile'ee 1 - If' + EL'^ ^e') vanishes. But then the hypothesis of Case 2 implies 

n r ° = n t ° = °- ( 7 - i? ) 

and again we conclude that at least two first order even theta constants vanish. This shows 
that the variety cut out by the polynomials r, s,,t, is contained in the Torelli boundary. 

We now show the reverse containment. Suppose a point p lies in the Torelli boundary. 
By [26, Theorem 3.1] at least 6 first order even theta constants T e e , vanish at p. In order to 
show that all 15 polynomials r, s,,t, vanish at p, we divide our analysis into three cases: 

Case 1. Suppose that at least two T ° , T° vanish. Then there exists e 6 Fj\ {0} such 

that one of e^e^ is in e and the other is not. For this choice of e, equation (7.17) holds. 
Since the expressions (7.10) and (7.11) for r and the s a each contain a factor of either 

(IL'Ge-L T e' + EL'^ T e') 01 (IL'Se-L T e> ~ Yle'fe^- T e') > We Se<3 tllat T aI1 d 3,11 S « Vanish. AlsO, 

each t a vanishes because the expression (7.12) has a factor ri e 'eF>j 

Case 2. Suppose that exactly one T® vanishes. Since at least six theta constants vanish, 
at least 5 of the T e e , with e ^ vanish. So there exist two vanishing T 6 , 1 , T e , 2 with ei, e 2 7^ 

and €j 7^ e 2 . By Riemann's theta relations (7.14), (Yl e i &e ±T®, — Yl^^T®} vanishes for 
e = ei,e2- Since the expression (7.10) for r contains all 7 such expressions as factors and 
each s a contains all but one of these as factors, we see that r and all the s a vanish. The t a 
vanish for the same reason as in Case 1. 

Case 3. Suppose that none of the T° vanishes. Then r and the s a vanish by the argument 
in Case 2. Using (7.16), we rewrite equation (7.12) for the t a as follows: 

U = C H n TyW a) (7.18) 

e£o- x e'Ge- 1 - 

where C is nonzero. Let T e /, j = 1, . . . , 6 be vanishing first order even theta constants whose 

characteristics form an azygetic 6-set. The vanishing of each of the TV implies the vanishing 

of 4 of the t a , namely those with a ^ ef. If the F2- vector space spanned by the €j is the 
whole space Fg, then all 7 t„ vanish in this way. If not, without loss of generality, we may 
assume that {e±, . . . , e 6 } C {001, 010, 011}. In this case, it suffices to show that ti o vanishes. 

We note that there are many ways to write W a as a polynomial in the T e e ,, coming from 
different ways of grouping the 4 theta constants in (7.9) into two pairs. For W wo , we have 

U/ — {tOOI \ 2 (rp001\2 /t-i001\2/t-i110\2 

Wioo — 000 ) \ 1 ioo) ~l J oioJ K 1 001 ) 

_ /t-tOIOn 2/^-1010^2 _ /t-i010\2/ati101\2 

— l J 000 J \ 1 W0) \ 1 001 I 

_ (rp01l\2(rp01l\2 _ ( rp01\\2 ( rp\l\\2 

— l J 000 / l J 100/ I J 011 / I J 011 / • 

By our assumption, six of the twelve T-variables in the expression for W100 vanish. The only 
way W\qq can be nonzero is for the following three conditions to hold: 

t-iOOI _ t^OOI _ n t^OOI _ rp\W _ n . 

J ooo ~~ J ioo ~~ u U1 J oio — J 001 ~~ u ) 
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rpOW 

1 000 

7-1011 
J 000 



mOlO 
1 100 

rpOU 

1 100 



or 

or T ° lT 



71101 

J 010 

rplll 

J 011 



This gives eight cases in total. It can be verified that in each of these eight cases, none of the 
characteristics of the six vanishing first order even theta constants are azygetic. We conclude 
that all t/j vanish, and this completes our proof that the Torelli boundary is contained in the 
common zero set of our coefficient polynomials r(u), So-(u), t a (u) for the Coble quartic. 

(c) Our next goal is to show that {F, S) is a prime ideal. Let C denote the subscheme of 
P 7 xP 7 defined by this ideal. We begin by proving that C is reduced and irreducible. Consider 
C as a family ir: C — > $(.4.3(2,4)). Let U denote the non-hyperelliptic locus in $(4.3(2,4)). 
The fiber over each closed point in hi is an irreducible Coble quartic hypersurface. Thus, 
7r _1 (W) is a family over an irreducible base whose fibers are irreducible and have the same 
dimension. Therefore, ir~ l (U) is irreducible by [19, Exercise 14.3]. Since F and S have no 
common factor, the ideal (F, S) is a complete intersection. Hence it is Cohen-Macaulay, and 
all of its minimal primes have the same dimension 12. 

We claim that C is the closure of tt (14) in P 7 x P 7 . Suppose it is not. Then it has a 
twelve-dimensional component contained in ($(4.3(2,4)) \U) X P 7 . Since dim($(4.3(2, 4)) \ 
U) = 5, this can only happen if there is a five-dimensional subvariety Z of $(4.3(2,4)) \ U 
such that Z x P 7 c C, i.e. all of the r,s a ,t a vanish on Z. This is impossible because the 
zero locus of r,s a ,t a is the Torelli boundary, which has dimension four. Therefore, C is 
irreducible and C = C. We now know that the radical of (F, S) is a prime ideal. 

The Cohen-Macaulay property implies that (F,S) satisfies Serre's criterion (Si). Since 
the variables x do not appear in S, the Jacobian matrix of (F, S) is a 2 x 16-matrix in block 
triangular form. If we pick any point (u, x) with u G S non-hyperelliptic and x non-singular 
on its Coble quartic, then the Jacobian matrix has rank 2 at this point. So, this point is 
non-singular on the scheme defined by F = S = 0, and hence (F,S) is generically reduced. 
Therefore, (F, S) is radical, and it is the prime ideal defining the universal Coble quartic C 
in P 7 x P 7 . The variety C has bidegree 16U(28U + AX) because it is a complete intersection 
defined by two polynomials of bidegree (16,0) and (28,4). This completes the proof. □ 

Remark 7.2. The Satake ring C[u]/ (S) is not a unique factorization domain (UFD). Our de- 
gree 28 polynomial r(u) has distinct factorizations. For any a e Ff, we have r = c a q a , where 

«.= n (n^-n^) ^ ^nfn^-n^) 

Both c a and q a are polynomials in the square of the and rLeF 3 ^e°- As before, we 
can rewrite them as polynomials in u. Therefore, as polynomials in u, we have r = c a q a 
mod (S). Hence c a is a factor of r in C[u]/(«S). It can also be verified that these c a are 
pairwise relatively prime in C[u], thus pairwise relatively prime in C[u]/ (S) since deg(co-) = 
12 < 16 = deg(5). This is impossible for a UFD. Compare this to results of Tsuyumine [51] 
and others on the UFD property for coordinate rings representing the moduli space 4.3. □ 

We close with a remark that highlights the utility of the formulas derived in this paper. 
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Remark 7.3. Our supplementary files enable the reader to write geometric properties of 
plane quartic curves explicitly in terms of the theta constants (3.4). For example, consider the 
condition that a ternary quartic is the sum of five reciprocals of linear forms. Such quartics 
are known as Liiroth quartics. A classical result of Morley states that Liiroth quartics form 
a hypersurface of degree 54 in the P 14 of all quartics. See [44] for a modern exposition. 
Liiroth quartics are also characterized by the vanishing of the following Morley invariant 

/l234567+/l234576 + /l234657 + /l234675 + /l234756 + /l234765 + /l235467 + /l235476+/l235647 + /l235674 
+/l235746+/l235764 + /l236457 + /l236475 + /l236547 + /l236574 + /l236745 + /l236754 + /l237456+/l237465 
+ /l237546+/l237564 + /l237645+/l237654 + /l243567+/l243576+/l243657 + /l243675 + /l243756+/l243765 

This expression is found in [44, page 379, after Figure 1]. Using the transformation derived 
in (6.11), the Morley invariant translates into the following linear form in Coble coefficients: 

6r + Sqoi + S 010 + s 100 + s 011 + s 101 + s 110 + s lll- 

We note that, for S 7 C W^(E 7 ) embedded as a parabolic subgroup, the 15- dimensional 
space of Coble coefficients decomposes as a 14-dimensional irreducible SV-module plus the 
trivial representation. The latter is spanned by the Morley invariant. Now, substituting the 
polynomials r(u) and Syfc(u) from Theorem 7.1 into this linear form, we obtain a polynomial 
of degree 28 that has 59256 terms. That expression in Mooo, Uooi, • • • , «iii represents the 
condition that a matrix r in the Siegel upper halfspace Sj^ comes from a Liiroth quartic. □ 

8 Equations for universal Kummer threefolds 

We now turn to the object that gave our paper its title. The Kummer threefold is the 
singular locus of the Coble quartic. In the past sections we found the ideals for three 
variants of the universal Coble quartic. Each is a twelve- dimensional projective variety, over 
a six-dimensional base. First, in Section 4, the base was P 6 . Next, in Corollary 5.3, the base 
was the Gopel variety Q . Finally, in Theorem 7.1, the base was the Satake hypersurface S. 

For each of the 3 versions of the universal Coble quartic we have a universal Kummer 
threefold. Each of them is a 9-dimensional irreducible variety. Their respective ambient 
spaces are P 6 x P 7 , P 14 x P 7 , and P 7 x P 7 . We discuss their defining equations in this order. 

First, consider the parametrization of the Coble quartic (2.12) in terms of ci,C2, ...,07 
given in (4.1). This defines what we call the flex version of the universal Kummer variety: 

^flcx c p6 x p7> 

For (ci : • • • : Cj) not lying on any of the reflection hyperp lanes of E7, we consider all points 
(c, x) = ((ci : • • • : C7), (xooo : ■ ■ • : £111)) such that x lies in the Kummer variety defined by 
c. The variety /C| ex is the Zariski closure of the set of all such points (c, x) in P 6 x P 7 . 

The label 'flex' refers to the fact that general points c in P 6 , modulo the action by W(Ej), 
represent plane quartic curves with a distinguished inflection point, as seen in the proof of 
Theorem 4.1. We shall construct two sets of polynomials for the ideal of /C| ex . The first set 
consists of eight polynomials. They are the partial derivatives, with respect to xooch • • • , %m, 
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of the Coble quartic (2.12). Here r, s,, t, are the degree 7 polynomials in ci, . . . , C7 listed after 
(4.1). Each of these Coble derivatives is the sum of 323 terms of bidegree (7,3) in (c,x). 

Our second set of equations for /Cg cx consists of 70 polynomials of bidegree (6, 4) in (c, x). 
Their construction is considerably more difficult, and we shall now explain it. The idea is to 
use the methods for degeneracy loci arising from Vinberg's ^-groups, as developed in [29]. 

We start with fi 3 (4), the third exterior power of the cotangent bundle of P 7 twisted by 
0(4). Its global sections are homogeneous differential 3-forms of degree 4 on P 7 . The space 
of global sections H°(P 7 , fi 3 (4)) is isomorphic to /\ 4 C 8 via the map 

A C 8 -+ H°(P 7 , fi 3 (4)) , a, A a, ^a k ^a ^ ^ ^ * f> A A f k \ ~ a ffi * ^ * f e \ 
I > J +aj(aai A dak A dag) — ai{daj A da^ A da{) 

Here the basis {a±, (12, ... , as} is denoted as in Remark 4.2. Let U±, U2, ■ ■ ■ , Us be the corre- 
sponding affine open charts on P 7 . We write U = Spec(C[zi, . . . , Zi-i, z i+ i, . . . , Zs\), where 
z s = a s /cii. The restriction f2 3 (4)|[/ i is generated by the 35 forms dzj A dzk A dzg where 
j,k,£ 7^ %. The fiber of the rank 35 bundle f2 3 (4) over any point x G U is identified 
with /\ 3 C 7 . The value of an element a G f2 3 (4) on that fiber can be written as a linear com- 
bination of the forms dzj A dz& A dzg. The 35 coefficients are homogeneous linear expressions 
in x, or non-homogeneous linear expressions in z. We regard these as the coordinates of a. 

Consider the action of GL^C) on /\ 3 C 7 . As explained in [29, §2.6], this action has unique 
orbits of codimensions 1, 4, and 7, respectively. We denote the closures of these orbits by 
Oi, 04,07, to indicate their codimension in /\ 3 C 7 . The orbit closure 0\ is a hypersurface 
of degree 7. Its defining polynomial / generates the ring of SL7(C)-invariant polynomial 
functions on /\ 3 C 7 . We constructed the invariant / using the method in [36, Remark 4.4]. 
Its expansion into monomials has 10680 terms. It is presented in our supplementary files. 
The ideal for O4 is generated by the 35 partial derivatives of /. These are polynomials of 
degree 6 in 35 unknowns, and we shall now introduce a certain specialization of these. 

Consider a section v = C\h\ + ■ ■ ■ + c 7 h 7 G /\ 4 C 8 where hi, . . . ,h 7 are the tensors in 
Remark 4.2, and (ci : • ■ • : c 7 ) G P 6 is generic. The section v is a linear combination of the 
56 forms ddj A daj A da& where each coefficient is a monomial of bidegree (1, 1) in (c, x). On 
Ui we set a, = 1 and aj = Zj for j ^ i. Then 21 of the 56 summands vanish since da; = and 
daj = dzj for j ^ i. What is left is a linear combination of the 35 forms dzj A dz\~l\dzg, where 
each coefficient is a monomial of bidegree (1, 1) or (1,0) in (c,z). These 35 coefficients are 
the coordinates of v. Shortly, we shall plug them into the derivatives of the polynomial /. 

The following degeneracy locus has codimension 4 in the 7-dimensional affine space Uf 

{xGf/j I v(x) G 4 }. (8.1) 

It was shown in [29, §6.2] that its closure in P 7 is precisely the Kummer threefold over c. 
Here the coordinates on P 7 need to be relabeled a\ = Xooch • • • > a 8 — x m as m Remark 4.2. 
For completeness, we mention that the corresponding degeneracy locus for O7 is the set of 
64 singular points of the Kummer threefold, but we won't need to use this. 

By plugging the 35 coordinates of v into the partial derivatives of /, we obtain 35 poly- 
nomials in (c,z). These polynomials are homogeneous of degree 6 in c, and they are non- 
homogeneous of degree 4 in z. For generic c, these equations define the affine Kummer 
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threefold in Ui = {cij ^ 0} ~ A 7 . If we homogenize these equations, then we obtain 35 
bihomogeneous polynomials of degree (6, 4) in (c, x). These all lie in the prime ideal of /C| cx . 

We now repeat this process for the seven other affine charts Uj. This leads to 8 ■ 35 = 280 
polynomials of bidegree (6,4) in (c, x), but only 120 of them are distinct. Some of these 
polynomials have 332 terms, while the others have 362. They span a 70-dimensional vector 
space over Q, and we select a basis for that space. We conjecture that this basis suffices: 

Conjecture 8.1. The prime ideal of the universal Kummer threefold in P 6 x P 7 is minimally 
generated by the 78 polynomials above, namely, 8 of bidegree (7,3) and 70 of bidegree (6,4). 

Second, let us consider the Gopel version of the universal Kummer variety: 

/c| opcl c g x p 7 c p 14 xp 7 . 

Here, the universal Coble quartic gives an equation of bidegree (1,4). Its eight partial 
derivatives are polynomials of bidegree (1, 3). These polynomials are not sufficient to generate 
the Kummer ideal, even over a general point in Q, because we need 70 quartics. The eight 
derivatives only give 64 equations of bidegree (1,4), so we need at least six more polynomials 
of bidegree (?, 4). We do not know how to produce these extra quartics. In other words, we 
do not know how to lift the degeneracy locus construction of (8.1) to the Gopel variety Q. 
In light of the beautiful combinatorics in Section 6, it is desirable to study this further. 

Last but not least, we return to the object that started this project. The theta version 
of our variety is the Zariski closure of the image of the universal Kummer map k in (3.8): 

£ 3 (2,4) C S x P 7 C P 7 x P 7 . (8.2) 

As before, we use the coordinates u = (uooo : ^ooi : • • ■ : ^111) on the first P 7 to parameterize 
the moduli of Kummer threefolds, and the coordinates x = (xooo : ^001 : ■ • • : £111) on the 
second copy of P 7 to parameterize points of a particular Kummer threefold. The theta version 
of the universal Kummer threefold has codimension five, and we have already constructed 
several polynomials in its defining bihomogeneous prime ideal X 3 . First, the ideal X 3 contains 
the bidegree (16, 0) polynomial of the Satake hypersurface S. Second, since the Kummer 
threefold is the singular locus of the Coble quartic hypersurface, one has the eight partial 
derivatives of F as in Theorem 7.1 with respect to x^. These have bidegree (28,3). 

Third and most important, there are additional generators of bidegree (16,4) in (u, x). 
These play the same role as the 70 equations of bidegree (6, 4) in (c, x) for the flex version 
in Conjecture 8.1. We will show that there are 882 such additional minimal generators. 

Lemma 8.2. There exists a polynomial f of bidegree (16,4) in the universal Kummer ideal 
I3, having the same form (2.12) as the Coble quartic, but now r,s,,t, are polynomials of 
degree 16 in u. These can be given explicitly by the formulas 

f — s ooi = S100 = s ioi — s no = sni = £010 = £110 = 
s 010 = (^000^011 + %01 u 010 — W100W111 — ^101 M llo)( M 000 M 011 + %01 u 010 + u 100 u lll + ^ioi' u iio) 
(WOOOWOII - MOOIMOIO + WlOOWlll - 1*1011*110 ) (W000«011 ~ ^001^010 - UxooUm + U 101 Ui 10 ) 
(«000«010«10lWlll - UooiUoiiUwoU lw ){u 00 oU 01Q u W ou lw - U 0Q1 U 0n Wioi Win) 
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Soil = -OoOOWoiO + "001^011 + MlOOMllO + u wl u ul )(u 00 oUo 10 + Mooi^on - u 100 u 110 - u wl u ul ) 
(«000«010 - "OOlWOH + ^100^110 - U W1 U U1 )(U 000 U 0W - UooiUon - Uioo^llO + MlOl^lll) 

(mooo«oh«ioo«hi - u 001 u ow u W iu lw )(u OQO u ou u wl u lw - n oin io^ioo^iii) 

^001 = -(WOOOWOII - MOOIMOIO ~ ^100^111 + Wl01«Ho) (^000^011 - MOOIMOIO + ^100^111 - WlOlWllo) 
(^000^011 + ^001^010 + UlOOWlll + mioi«ho)(«oooMoii + ^ooi^oio - uioo«m - «ioi«iio) 

/ \/22 22 22 22\ 

(%00%10 M 101 M 111 — u 001 u 011 M 100 w 110j(. w 000 U 010 ~~ W 001 W 011 + M 100 M 110 — w 101 W lllJ 

toil = (WOOOWOIO - MOOIMOU - MiooMllO + Wl0lWlll)(w000W010 - MOOIMOU + MlOOMllO ~ Wl01«lll) 

(^ooo^oio + ^ooi^on - ^100^110 - wioiWih)(mooo«oio + ^001^011 + MlOOMllO + Ul01«lll 

(«ooo«ou«ioi«uo - wooi«oio«ioo«m) ( U 000 U 011 ~ U 001 n 010 + U 100 U 111 ~ u ioi n iio) 

hoi = -(moooMou + MooiMoio + uioowm + mioiMiio)(moooMoii - MooiMoio + WiooWm - WioiWuo) 

(^ooo^on + ^ooi^oio - Wioo«iii - mioi "no) (^000^011 - ^ooi^oio - u 100 u in + Uioiiino) 

(n oo M oio u ioo u iio — M ooi n oii u ioi u in)(' u ooo' u oio ~~ n ooi n on ~~ ^ioo^iio + u ioi u m) 

hll = (WOOOWOIO - MOOIMOU + MlOO^HO ~ ^lOl^m) ("000^010 - MOOIMOU ~ Wl00«110 + WlOl^lll) 
OoOO^OlO + MoOlMoil + MlOOMllO + U W1 U U1 ) (u 000 U 010 + WoOl^OH ~ ^100^110 ~ 1*101^111 ) 
(moooMohMiooMih - M 0lM010Ml0lMll0)(«000 w 011 ~ M 001 M 010 ~ M 100 M 111 + w ioi u no) 

. _ r 4 4 4 4 "I _. r 6 2 4 2 2 ] , qT 4 2 2 4 2 2] 

r 100 — P000' u 010' u 100' u lllJ4 i L w OOO u 010"011 W 100"lloJl6 ' ^PoOO"001"010"011"100"lllJl2 

[4 4 4 4 i _ r 6 4 2 2 2] _ Q[ 4 2 4 2 2 2] 

[u 000 u 011 u 100 ti 110 j4 [ u ooo u oio M oii w ioo w mJi6 z L u ooo n ooi n oio n on M ioo M noJi2 
+4[u ^uooin \ u 01 iu 10 ou 101 ui 10 u m ] 4 - 4[no 00 uooin io^oii u ioo^ioi^iio'«iii]4 

The last coefficient tioo(u) is irreducible and has 72 terms. In particular, f is the sum of 
1168 monomials of degree (16,4) in (u, x). 

Bert van Geemen informed us that the polynomial / above admits a determinantal rep- 
resentation similar to (1.1). That determinant will be derived in a forthcoming paper of his. 

Proof. We start with the following quartic relation among the theta constants of genus four. 

0[OOlO|OOOl](r; 0) • 0[OO1O|1OO1](t; 0) • 0[OO1O|O1O1](t; 0) ■ 0[OOlO|llOl](r; 0) 

- 0[OO11|OO11](t; 0) ■ 0[OOll|lOll](r; 0) • 0[OOll|Olll](r; 0) ■ 0[OOll|llll](r; 0) 

- 0[OOlO|OOOO](r; 0) • 0[OOlO|lOOO](r; 0) • 0[OOlO|OlOO](r; 0) ■ 0[OOlO|llOO](r; 0) 

+ 0[OOll|OOOO](r; 0) ■ 0[OOll|lOOO](r; 0) • 0[OOll|OlOO](r; 0) ■ 0[OOll|llOO](r; 0) = 0. 

Genus four identities of this form can be derived from genus two theta relations using [24, 
Proposition 4.18]. Such relations hold identically in r G i^4, not just on the Schottky locus. 

The main computation is to express the relation (8.3) in terms of the genus four moduli 
variables u. To do this, we first turn (8.3) into a relation between squares of the 0[e|e'](r; 0). 
This can be done by taking the product of the relation (8.3) with its seven conjugates, where 
all possible signs in front of each term are chosen. This gives a polynomial of degree 16 in 
the 2 [e|e'](r; 0). We next apply the genus four version of (3.6). This writes the squares of 
the first order theta constants as quadrics in the 16 moduli variables Moooo5 M oooi ? • • • j w mi- 
The result is a huge homogeneous polynomial of degree 32 in these 16 variables. 

To this polynomial we now apply the Fourier-Jacobi expansion technique as in (3.5). 
This replaces the polynomial above by one of its initial forms, but now in the 16 unknowns 
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Uijk,Xijk- The result is a polynomial of degree (28,4) having 12268 terms. It lies in the 
universal Kummer ideal X 3 and is expressible as a Q[u] -linear combination of the 15 degree 
four invariants. Remarkably, this degree (28, 4) polynomial turns out to be reducible. It has 
an extraneous factor of degree (12, 0). That factor is a polynomial in u which cannot vanish 
on S. Dividing out this factor gives the desired bidegree (16,4) polynomial / £ X 3 . □ 

The action by the modular group Sp 6 (Z) on C[u, x] induced by the action on P 7 x P 7 
preserves the ideal X 3 . On the space of Heisenberg invariant polynomials, the principal con- 
gruence subgroup T 3 (2) acts trivially, so there is a well-defined action of the quotient group 
Sp 6 (F 2 ) = Sp 6 (Z)/r 3 (2), cf. (3.1). For a detailed discussion of the relevant representation 
theory, see [14]. However, the polynomial / in Lemma 8.2 is not invariant under this action. 
Therefore, we obtain more polynomials in X 3 by applying this action to /. 

Lemma 8.3. The orbit of the polynomial f in Lemma 8.2 under the action of Sp 6 (F2) 
contains exactly 945 elements. They span a vector space over C of dimension 882. 



Proof. In what follows, we will use a right action of Sp 6 (F2). Recall that on the projective 
coordinates (u, x), generators of Sp 6 (F 2 ) of the form 71(A) = (A' 1 )*^ aC ^ ^ 

71(A) o a a = a^o for a £ {u,x}. 

Generators 72CB) = with B f = B act by 

/ \ 7TZ<T tig r ~. 

7 2 (-B) o a a = e 2 a a for a £ [u,x\, 
and the Weyl element 7 3 = ^ ^\ acts by the finite Fourier transform 

73°o CT = ^(-l) p<T a p for ae{u,x}. 
Alternatively, the group Sp 6 (F 2 ) is generated by the following two elements // and v'\ 

(\ o\ (\ 0\ 

// = j 2 {B) 7l {A) for B = , A = 1 , (8.4) 




v' = 71(^)72(^)7372(5)7372(5) for A = 1 0,5= . (8.5) 



These two generators correspond to \x and v in (4.2). They act on C[u, x] by the substitutions 

/ a ooo l— >* a oocb a ooi ^ °ooi) °ioo i ' a ioi 5 a ioi ^ i ' ^ioo; t ^ r • / — t 

jj, : . tor a £ {u,xf, % = y— 1, 

aoio l— >■ ooiO) oon h ~>* a oii ; a no ^ * ■ «m, iui l— >■ * ■ ono 

/ aooo ^ flooo+flioo; flooi ^ flooo - flioo; floio ^ a ooi+flioi, O011 1— )■ aooi— aioi t _ r n 
z/ : tor a £ {w, x}. 

aioo | — ► Ooio+flno, 0101 a oio — Oho, Quo ^ %1+ain, am H- aon - om 
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We can try to generate the orbit by applying the generators // and v' successively. This 
is challenging because some of the polynomials in the orbit are very large. Instead, for some 
N > 2, we choose random vectors (uj,Xj) over a finite field, and we evaluate 

(/(a(ui,xi)) : ■ ■ • : /(a(ujv,Xiv))), (8.6) 

where a is the product of a sequence of /i's and v's. In this way, we get 945 distinct points 
in P^ -1 . Therefore, the orbit of / under Sp 6 (F2) has at least 945 elements. 

To prove that the number 945 is exact, we study the stabilizer of / in Sp 6 (F2). It suffices 
to show that the stabilizer has at least |Sp 6 (F 2 )|/945 = 1536 elements. It can be seen from 
the explicit form of / in Lemma 8.2 that the stabilizer contains 71(A) for the 8 unipotent 
lower triangular 3 x 3-matrices A, and all of the 64 elements 72 These elements generate 
a group of order 8 • 64 = 512. Note that all of these elements act on P 7 x P 7 by signed 
permutations on the coordinates. We found another element in the stabilizer, namely 



1/71 (A) for A= 1 . (8.7) 




This element of order 3 does not act by signed permutations on the coordinates and thus is 
not in our subgroup of order 512. Since the size of the stabilizer divides |Sp 6 (F 2 )| = 1451520, 
it must be at least 512-3 = 1536. This shows the orbit has size 945. Alternatively, we can do 
explicit calculations with these matrices in GAP and establish the exact size of the subgroup. 

A linear algebra computation was used to show that the vector space spanned by the 
orbit of / has dimension > 882. We generated the 945 polynomials of degree (16,4) in 
implicit form as in (8.6). We evaluated them at random configurations of N > 945 points 
over various finite fields. In each experiment, the resulting 945 x iV- matrix has rank 882. 
Hence the C-vector space spanned by the 945 orbit elements has dimension at least 882. 

To prove that the number 882 is exact, we identified exactly 945 — 882 = 63 linearly 
independent relations among the 945 orbit elements. Fortunately, the relations are rather 
simple. There exist 15 distinct group elements gi, ■ ■ ■ ,gi5 G Sp 6 (F 2 ) such that 

15 

5>/°» = °- (8-8) 
i=i 

The complete list is found in our supplementary materials. Applying Sp 6 (F 2 ) to (8.8) gives 
63 linearly independent relations since no two relations involve the same conjugate of /. □ 

Remark 8.4. We can verify that our polynomials lie in X 3 using direct numerical compu- 
tations. Indeed, by running Swierczewski's code for theta functions in Sage [50], we can 
generate arbitrarily many points (u, x) on the universal Kummer variety £3(2,4). The 
polynomial with 1168 terms was shown to vanish on all of them. Likewise, we can check 
numerically that this vanishing property is preserved under the two substitutions //, v' . 

Remark 8.5. The space spanned by the Sp 6 (F 2 )-orbit of / is a representation of Sp 6 (F 2 ). 
It is equivalent to a subrepresentation of the representation induced from the trivial repre- 
sentation of the stabilizer subgroup of the polynomial / from Lemma 8.2. According to a 
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GAP calculation, the irreducible subrepresentations of that 945-dimensional induced Sp 6 (F 2 )- 
representation have the following dimensions: 

1, 27, 27, 35, 35, 84, 120, 168, 168, 280. 

Here the two occurrences of 27 and 35 refer to an irreducible representation that appears 
with multiplicity 2. Note that 945 — 882 = 63 = 1 + 27 + 35, and this is the unique way to 
build 63 from the dimensions of these irreducible summands. So we can use this to determine 
the structure of these 882 equations as a representation of Sp 6 (F 2 ). □ 

We do not completely understand what happens to the universal Kummer threefold 
£3(2, 4) and the generators found above when we restrict the base to one of the 36 hyperel- 
liptic divisors in S. The issue is that, for fixed r in the hyperelliptic locus, the ideal of the 
Kummer threefold has a different structure. It can be seen from equations (7.10), (7.11) and 
(7.12) that the vanishing of a theta constant T e e , causes the universal Coble quartic to become 
the product of a monomial in the remaining non-zero theta constants times the square of a 
quadric with integer coefficients in the x variables. For instance, if Tq = then 

F T = (^000 + 3^001 + X 010 + X 011 + X 100 + X 101 + X 110 + X lll) ' Y\. C^e') ■ (8-9) 

eG Fl\{0} 

Hence over the hyperelliptic locus, the prime ideal of a Kummer threefold contains a quadric. 
This can also be seen from the vanishing of a theta constant 0[e|e'](r;O) in equation (3.9). 
That quadric generates 8 linearly independent cubics and 36 linearly independent quartics. 
The additional 34 quartic relations were identified recently by Miiller [41, Theorem 4.2]. 

We optimistically conjecture that the equations we have found so far are sufficient: 

Conjecture 8.6. The prime idealT^ of the universal Kummer variety (8.2) is generated by 
891 bihomogeneous polynomials in (u, x); the Satake polynomial of degree (16,0) , the eight 
Coble derivatives of degree (28,3), and the 882 polynomials of degree (16,4). 

What remains to be done is to better integrate the generic case and the hyperelliptic 
case. This will be crucial for understanding the relationship between our 70 extra quartics 
over P 6 and our 882 extra quartics over S. How can we construct these over Ql How can 
we lift the map (7.1) to the universal Kummer threefolds? We expect that the study of 
hyperelliptic moduli in [21] and its Macdonald representation of S$ will be relevant here. 

In closing we remark that the Fourier-Jacobi method was the key to success in Lemma 8.2. 
Here the systematic passage to a non-trivial initial form was driven by a toric (or sagbi) 
degeneration as in (3.5). This is fundamental also for tropical geometry, our next topic. 

9 Next steps in tropical geometry 

We now take a look at Kummer threefolds and their moduli through the lens of tropical 
algebraic geometry [9, 10, 16, 30, 31, 40, 48]. Each of our ideals defines a tropical variety, 
which is a balanced polyhedral fan. These fans represent compactifications of our varieties 
and moduli spaces, and they allow us to understand what happens when the field C of 
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complex numbers gets replaced by a field K with a non-trivial non-archimedean valuation. 
This section serves as a manifesto in favor of explicit polynomial equations. They are essential 
for understanding the combinatorics that links classical and tropical moduli spaces of curves. 

We begin our discussion with the six-dimensional Gopel variety Q. Recall from The- 
orem 5.1 that Q sits inside P 14 where it is defined by 35 cubics and 35 quartics. The 
tropicalization of this irreducible variety is a pure six- dimensional polyhedral fan in TP 14 = 
R 15 /M(l, 1, . . . , 1). This fan can be computed, at least in principle, with Gf an [34]. However, 
it does not have good combinatorial properties. For tropical geometers, it is much better 
to pass to the modification arising from the re-embedding, seen in Section 6, of the Gopel 
variety Q in P 134 . Its ideal is generated by binomials and linear trinomials. We define the 
tropical Gopel variety to be the tropical variety of that ideal. This is a six-dimensional fan 
which lives in TP 134 . The Weyl group W(E7) acts on trop(£?) by permuting coordinates. 

We shall construct the tropical Gopel variety trop(^) combinatorially, not from its defin- 
ing ideal via Gfan [34], but directly from the parametrization given by (6.14) and (6.15). 
This is equivalent to the parametrization in Section 4, but we have factored it as follows: 



The first map £ is given by evaluating the 63 linear forms (6.14) that represent E 7 . The second 
map m is the monomial map (6.15). It corresponds to the 63 x 135-matrix A that encodes 
incidences of vectors and Lagrangians in (F2) 6 \{0}. The tropicalization of the monomial 
map m is given by the classically-linear map TP 62 ->■ TP 134 defined by its transpose A 1 . The 
tropicalization of the linear space image(£) is the Bergman fan of the matroid M(E 7 ) of the 
reflection arrangement of type E 7 , as defined in [2]. This fan lives in TP 62 , and we have: 

Proposition 9.1. The tropical Gopel variety trop(£?) coincides with the image of the Bergman 
fan of M (E 7 ) under the linear map TP 62 — > TP 134 given by the matrix A 1 from Section 6. 

Proof. The factorization (9.1) shows that Q is the image of a map whose coordinates are 
products of linear forms. The result then follows immediately from [16, Theorem 3.1]. □ 

The first step towards the tropical Gopel variety trop((?) is to list the flats of the matroid 
M(E 7 ). In Table 2, we present the classification of all proper flats of the matroid M(E 7 ). 
There are 30 orbits of flats under the W / (E 7 )-action. For each iy(E 7 )-orbit we list the rank, 
the size of the orbit, and a set of linear forms that serves as a representative. These linear 
forms define an intersection of the hyperplanes in the reflection arrangement, of codimension 
equal to the rank, and the flat consists of all hyperplanes that contain that linear space. 

In Table 2, we are using the bijection between flats and parabolic subgroups of the Weyl 
group, which can be found in [3, Theorem 3.1]. These parabolic subgroups correspond to 
root subsystems in the finite type Dynkin diagram. We calculated Table 2 from scratch. 
Similar information for the lattice of parabolic subgroups can be found in [23, Table A. 2]. 
This can also be taken as an independent verification of the correctness of Table 2. 

Eleven of the 30 orbits consist of irreducible root subsystems. Their numbers are marked 
in bold face. The total number of proper irreducible flats of the matroid M(E 7 ) is therefore 




(9.1) 



6091 



63 + 336 + 1260 + 2016 + 315 + 1008 + 336 + 378 + 288 + 63 + 28. 
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Family # 


Rank 


Size 


Root subsystem 


Equations of a representative flat 


1 


1 


63 


Ax 




2 


2 


336 


A 2 


Or, C\ + c 3 — cq 


3 


2 


945 


Ai x Ai 


C 4 , C 7 


4 


3 


1260 


A 3 


c 7 , c 4 , ci - c 5 


5 


3 


5040 


Ai x A 2 


or, c 4 , c 3 + c 5 - c 6 


6 


3 


3780 


Af 


Oo, c 3 - c 5 , c 2 - c 7 


7 


3 


315 


Af 


or, c 4 , c 3 


8 


4 


2016 


A 4 


c 7; c 4i C 3 + C5 — C6, C\ — C 5 


9 


4 


315 


D 4 


c 7; c 5; c 4> Ci 


10 


4 


7560 


Ai x A 3 


c 7; c 5; c 4; c 3 — c 6 


11 


4 


3360 


A 2 x A 2 


C7, C 6 , C 3 - C 4 + C 5 , Ci + C 4 + C 5 


12 


4 


1260 


Ai x A 3 


c 7 , c 5 - c 6 , c 4 , c 3 


13 


4 


15120 


Af x A 2 


c 7) c 6) c 2 + c 3 + 2c 4 + C5, Ci + C 4 + C5 


14 


4 


3780 


Af 


c 7) c 6; c 5; c 4 


15 


5 


336 


A 5 


c 7) c 6) c 5; c 2 — C 4 ; Ci + C 4 


16 


5 


1008 


A 5 


c 7; c 6) c 4; c 3 + c 5) c l + c 5 


17 


5 


378 


D 5 


C 7 , C 6 , c 3 , c 2 - C 4 , Ci 


18 


5 


6048 


Ai x A 4 


C 7 , C 6 , c 4 , c 2 - C 3 - C 5 , Ci + C 5 


19 


5 


945 


Ai x D 4 


C 7 , C 6 , C 5 , C 4 , C 3 


20 


5 


5040 


A 2 x A 3 


c 7; c 6; c 4; c 2 — 2Cs, Ci + C5 


21 


5 


7560 


Af x A 3 


c l, Cq, C5, C 4 , Ci + C 2 


22 


5 


10080 


A 1 x Af 


c 7i c 6; c 4; C 2 — C 3 — C5, Ci — C 3 — 2Cs 


23 


5 


5040 


Af x A 2 


C 7 , Cq, C 4 , C 3 , Ci + C 2 - C 5 


24 


6 


288 


A 6 


C 7 , C 6 , C 4 , C 3 + c 5 , c 2 - 2c 5 , Ci + C 5 


25 


6 


63 


D 6 


C 7 , C 6 , C 5 , C 4 , C 3 , C 2 


26 


6 


28 


E 6 


C7,C 6 ,C 4 - C 5 ,C 3 ,C 2 - C 5 ,Ci 


27 


6 


1008 


Ai x A 5 


c 7; c 6; c 5; c 4; c 2 + C 3 , Ci — C 3 


28 


6 


378 


Ai x D 5 


c 7; c 6; c 5; c 4) c 2; Ci + C 3 


29 


6 


2016 


A 2 x A 4 


c 7; c 6; c 4; c 3 + c 5; c 2 — 2Cs, Ci — 3Cs 


30 


6 


5040 


Ai x A 2 x A 3 


or, c 6 , c 4 , c 3 , c 2 - 2c 5 , ci + c 5 



Table 2: The flats of the E 7 reflection arrangement. 
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The Bergman fan of M(E^) is a six-dimensional fan in TP . It has fi = 6091 rays, one for 
each irreducible flat F. The rays are generated by their 0-1-incidence vectors ep = ^2 ie p &i- 
A collection J 7 of flats is nested if, for every antichain {F\,F2, . . . ,F r } in J 7 with r > 2, 
the flat F\ V F2 V ■ ■ ■ V F r is not irreducible. (This flat represents the subspace obtained by 
intersecting the given subspaces). For any nested set J 7 , we consider the convex cone Cj 
spanned by the linearly independent vectors ep, F G J 7 . The Bergman fan of M(E 7 ) is the 
collection of all cones Cjr where J 7 runs over all nested sets of irreducible flats of M(E 7 ). 

By the results of Ardila, Reiner and Williams in [2, §7], this simplicial fan is the coarsest 
fan structure on its support. Further, if fi is the number of i-dimensional cones Cj- then 

h - h + h - h + /a - fx + 1 = 1 ■ 5 ■ 7 ■ 9 ■ 11 ■ 13 • 17 = 765765. (9.2) 

Equation (9.2) rests on two non-trivial facts from matroid theory. First, the simplicial 
complex underlying the Bergman fan is a wedge of fi spheres, where fi is the Mobius number 
of the matroid (see [1, p. 42, Corollary]). Second, for the matroid of a finite Coxeter system 
(as in [2]), the Mobius number fi is the product of the exponents of that group [43, (1.1)]. 

We note that Bramble [7] was the first to determine the fundamental invariants of W(E 7 ). 
As for any reflection group, their degrees 2, 6, 8, 10, 12, 14, 18 are the exponents plus one. 

Corollary 9.2. The tropical Gopel variety trop(^) in T 134 is the union of the convex poly- 
hedral cones A 1 Cjr where J 7 runs over all nested sets of irreducible root subsystems 0/E7. 

Proof. This follows from Proposition 9.1 and the construction of the Bergman fan in [2]. □ 

Another approach to the tropical Gopel variety trop(^) is to use the formulas (6.6) and 
(6.8) for the Gopel functions f 9 ,g, as polynomials of degree 7 in the brackets [ijk]. This 
defines a rational map from the Grassmannian Gr(3, 7) C P 34 to the Gopel variety Q C P 134 . 
The tropicalization of this map is a piecewise linear map fi: TP 34 — > TP 134 . Note that the 
map fi is not linear because the expression for g, in the brackets [ijk] is a binomial and not 
a monomial. We expect that this rational map does not commute with tropicalization. 

The image of the tropical Grassmannian trop(Gr(3, 7)) under the piecewise linear map fi 
is a subfan of the tropical Gopel variety trop(^). It would be very interesting to identify that 
subfan. Note that, by [31, Theorem 2.1], that tropical Grassmannian has the face numbers 

/(trop(Gr(3,7))) = (721,16800,124180,386155,522585,252000). 

The tropical Gopel variety trop(^) is a modification and its face numbers are even larger. 

The positive part of the tropical Grassmannian governs the combinatorics of the cluster 
algebra structure on the coordinate ring of Gr(3, 7). Interestingly, the Weyl group relevant 
for this is of type E$ and not E7. Namely, it is shown in [48, §7] that the normal fan of 
the E 6 -associahedron defines a simplicial fan structure on the positive part of trop(Gr(3, 7)), 
with f- vector (42, 399, 1547, 2856, 2499, 833). The relationship to cluster algebras is explained 
in [48, §8]. In light of this, it would be interesting to examine the positive part of trop(^). 

Here is another question: What is the precise relationship between the Gopel variety Q, 
or the Satake hypersurface S, with the compactified moduli space Y(A) of degree 2 del Pezzo 
surfaces that was constructed by Hacking, Keel, and Tevelev in [30]? The latter is a tropical 
compactification of the space of configurations of seven points in P 2 . The underlying tropical 



14 



variety is a six- dimensional simplicial fan .F(A) that lives in an irreducible H A (E7)-module of 
dimension 35 [30, Corollary 5.3]. According to [30, §1.14], the link of J 7 (A) is the simplicial 
complex of nested sets on the root subsystems having types A 1; A 2 , A 3 x A 3 , and A 7 . 

The Satake hypersurface S lives in P 7 , and it is parametrized by fjs via the theta constant 
map The Newton polytope of its defining polynomial of degree 16 has face numbers 

/ (Newton (S)) = (344,2016,3584,2828,1120,224,22). 

The corresponding tropical hypersurface trop(iS) in TP 7 is a six-dimensional fan with 22 rays 
and 2016 maximal cones. Now we shall explain and derive the following result: 

Proposition 9.3. The image of the tropical Siegel space trop($) 3 ) under the piecewise-linear 
map trop(?9) is the intersection o/trop(iS) with the normal cone 0/ Newton(5) at the vertex 

M = -2«o 00 uooiWoio«ohWioo«ioiMiio«iii- (9-3) 

The map trop(^) induces the level structure on trop(^4.3) described by Chan in [9, §7.1/. 

Proof. We first define the terms in Proposition 9.3. Following [10, 40], the tropical Siegel 
space trop(i^3) is the cone PD3 of positive definite real symmetric 3 x 3-matrices. We use 
tropical theta functions as described by Mikhalkin and Zharkov in [40]. The coordinates of 
the tropical theta constant map trop(?9) are indexed by a £ {0, l} 3 , and they are defined by 

trop(tf) CT (T) = mm{(n + cr/2)* ■ T ■ (n + a/2) : n £ Z 3 } for T £ PD 3 . (9.4) 

The function trop(-$) CT from PD 3 to M is well-defined and takes non- negative values since T 
is positive definite. It is zero for all T when a = (0, 0, 0). Hence we have a well-defined map 

trop(tf) : PD 3 trop(S) C TP 7 . 

To show that the image lands in the tropical Satake hypersurface, we set r = ipT in (3.3) 
and (3.4), where p — > 00 is a real parameter. Write 6 2 [cr](r; 0) as a series in e = exp(— p), 
with (9.4) as the exponent of the lowest term. Then each coordinate of i?(t) is a series in e. 
Plugging these eight series into the Satake polynomial S, we obtain zero. Compare the 471 
monomials of S according to their order in e after this substitution. Two (or more) of the 
monomials must have the same lowest order in e. This means that trop($(T)) lies in trop(«S). 

The vertex of Newton(iS) given by the monomial M in (9.3) is simple, i.e. it has precisely 
seven adjacent edges. These edges are the Newton segments of the seven binomials 

M 000 U 001 W 010 W 100 M 111 + ^OOO^OOI^OIO^IOI^IIO + Wooo M 001 M 011 M 100 M 110 + 

u ooo tt ooi' u oii' (i ioi' u iii + Mooo^oio^on^ioo^ioi + w 000 m 010 m 011 m 110 m 111 + M, (9.5) 

U 000 U 100 M 10l' lt 110 M lll + M. 

Our computations revealed that these seven are all the e-leading forms of S selected by generic 
matrices T £ PD 3 . Hence the image of trop(^) consists of the seven outer normal cones at 
these edges. This is precisely what we had claimed in the first statement in Proposition 9.3. 

The correspondence with the level structure described in [9, §7.1] is seen as follows. 
The domains of linearity of trop($) define a subdivision of PD3 into infinitely many convex 
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Codim 


Glass# 


Geometric description 


Components 


How Many 


1 


1 


Hyperelliptic locus 


quadric HiS 


36 


2 


6 


Torelli boundary 


p3 x pi 


336 


3 


9 


Product or three elliptic curves 


P 1 X P 1 X P 1 


1120 


3 


16 


Satake boundary *4. 2 (2,4) 




126 


4 


18 


Torelli boundary of A 2 (2,4) 


quadric in P 3 


1260 


5 


24 


Satake stratum Al(2,4) 


pi 


1260 


6 


28 


Satake stratum Aq(2,4) 


point 


1080 



Table 3: Glass strata in the Satake hypersurface 



polyhedral cones. This subdivision is a coarsening of the second Voronoi decomposition. We 
group the cones into seven classes, according to which binomial in (9.5) gets selected by T. 

The seven classes in PD3 are naturally labeled by the seven lines in the Fano plane P 2 (F2). 
These lines are given by the three variables missing in the respective leading monomials in 
(9.5). For instance, the first binomial in (9.5) determines the line {(1:1:0), (1:0:1), (0:1:1)} 
because wno, 1*101, ^on are missing in the monomial prior to M. This subdivision of PD 3 into 
seven classes is precisely the level structure that was discovered by Chan in [9, §7.1]. □ 

We expect an even more beautiful structure when tropicalizing the re-embedding S' of 
the Satake hypersurface S into P 35 given by the quadrics (T^,) 2 in (7.13). To see this, we 
revisit the combinatorial result in [26, Theorem 3.1]. Glass classified the points in S' C P 35 
according to which of the 36 coordinates are zero at that point. Up to symmetry, there are 
seven Glass strata in S'. Table 3 identifies the irreducible components of these strata and 
their geometric meaning. Glass# is the number of coordinates (T e e ,) 2 that are zero. 

Table 3 refines our analysis of the Satake hypersurface in Section 3. For each irreducible 
component we know their defining polynomials in the unknowns uqqq, uqqi, . . . , «m. These 
played an important role in this paper. For instance, the components for Glass# 16 are the 
spaces Hf ~ P 3 in (7.2), and the components for Glass# 6 are the Segre varieties in (3.7). 

Now consider the tropical variety trop(iS') in TP 35 . This is a modification of the tropical 
hypersurface in TP 7 given by Newton(iS). It has much better properties since the coordinate 
functions now have a geometric meaning: they are the 36 divisors of the hyperelliptic locus. 
These appear as distinguished rays in trop(iS'). By the principle of geometric tropicalization 
[30, §2], they span a distinguished subfan of trop(iS'). Here, a collection of coordinate rays 
spans a cone if and only if the corresponding intersection of divisors appears in Table 3. 

The poset of Glass strata seems to embed naturally into the poset of cones in the tropical 
moduli space trop(-4 3 ). Compare [26, Theorem 3.1] with [9, Figure 8]. This deserves to be 
studied in more detail. Can the table above be lifted to a tropical level structure on A3I The 
finite groups of Section 3 act on the fans we described. Taking quotients by these groups 
leads to objects known as stacky fans. The tropical Torelli map of [9, 10] is a morphism of 
stacky fans from trop(Ais) onto trop(^4a). It would be desirable to express the combinatorics 
of the tropical Torelli map directly in terms of the polynomials and ideals in this paper. A 
related question: if we lift the map (7.1) to S' C P 34 , how does its tropicalization look like? 
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